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Abstract 

We study a dynamical system consisting of a massive piston in a cubical con- 
tainer of large size L filled with an ideal gas. The piston has mass M L!^ 
and undergoes elastic collisions with N ^ non-interacting gas particles of mass 
m = 1. We find that, under suitable initial conditions, there is, in the limit L — > oo, 
a scaling regime with time and space scaled by L, in which the motion of the piston 
and the one particle distribution of the gas satisfy autonomous coupled equations 
(hydro dynamical equations), so that the mechanical trajectory of the piston con- 
verges, in probability, to the solution of the hydrodynamical equations for a certain 
period of time. We also discuss heuristically the dynamics of the system on longer 
intervals of time. 

Contents 

1 Introduction 

2 Hydrodynamical equations 

3 Dynamics before the first recollision 

4 Dynamics between the first and second recolhsions 

5 Beyond the second recohision 
Appendix 

BibHography 

""^Department of Mathematics, University of Alabama at Birmingham, Alabama 35294 
^Department of Mathematics, Rutgers University, New Jersey 08854 
•^Department of Mathematics, Princeton University, New Jersey 08544 
''Current address: Institute for Advanced Study, Princeton, NJ 08540 



1 



1 Introduction 



The evolution of a macroscopic system consisting of a gas in a container divided by a 
massive movable wall (piston) is an old problem in statistical physics with a colorful his- 
tory. It was discussed by Landau and Lifshitz [LL] and later by Lebowitz [LI], Feynman 
[F], Kubo [Ku], see recent surveys by Lieb [Li], Gruber [G] and others [KBM]. 




Figure 1: Piston in a cylinder filled with gas. 

In its simplest form, the model consists of an isolated cylinder filled with gas and 
divided into two compartments by a large piston which is free to move along the axis of 
the cylinder, see Fig. 1. Initially, the piston is held fixed by a clamp and the gas in each 
compartment evolves independently and is in equilibrium^. We denote the density and 
temperature in the left and right compartments by tii^Tl and nji,Tji, respectively. The 
gas exerts pressure (= force per unit area) on the piston, which is given by equilibrium 
statistical mechanics as a function of density and temperature: Pl = PinL^Ti) and 
Pfj = PijiRiTR) on the left and on the right, respectively. At time t = the clamp 
is removed and the piston is released. Now one wants to describe the evolution of the 
system, especially its limit (final) state as t — > cxd. 

Starting with Pl 7^ Pr, the piston moves under the net pressure difference and 
compresses the gas whose pressure is lower, until its pressure builds up and it pushes 
the piston back. Depending on the initial values of n^, Tj^, ra^, Tr and the dynamical 
characteristics of the gases, the piston may follow a complicated trajectory, sloshing back 
and forth, but gradually it comes to rest at a place where the pressures are equalized on 
both sides: Pl = Pr- At that time one also expects that the gas in each compartment 
will again be in equilibrium. 

^For an isolated system of N atoms and total energy the equilibrium distribution is defined in 
statistical mechanics by a uniform probability distribution pcq on the energy surface E — const in the 
phase space (peq is called a microcanonical ensemble, it remains invariant under the dynamics by the 
Liouville theorem). One says that the system (gas) is "in equilibrium" if its states are "typical" for the 
measure pcq- In this state, a macroscopic gas will have (approximately, for large N) a uniform spatial 
density and a Maxwellian velocity distribution. The latter is defined so that the x, y, z components of the 
velocity vectors are independent normal random variables Af(0, ti^) with the same variance — ksT Im, 
where fc^ is Boltzmann's constant, T the temperature of the gas (which is a function of E, see below), 
and m the mass of an atom. 
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We observe, however, that the equahty of pressures Pl — Pr and the fact that the 
gas in each compartment separately is at equihbrium does not guarantee that — Tji. 
In particular, for dilute gases (which we shall consider from now on), the pressure is 
related to the density and temperature by P = nksT, where ks is Boltzmann's constant, 
and it is possible that Tl < Tr while ul > ur, so that the gas in the left compartment 
is cooler but denser, and in the right one hotter but more dilute (or vice versa). The 
exact values of the temperatures and T/j, at the time when there is no longer any 
pressure difference between the left and right and the piston comes to rest, depend on 
the initial conditions and other characteristics of the system, see [CPS, G], for example. 
Therefore, we have two possibihties now. If it happens that Tl = Tr as well as Pl = Pr, 
then the system as a whole will be in equilibrium, and we say that it came to a thermal 
equilibrium. On the contrary, if Pl — Pr but Tl ^ Tr, the system is said to be at 
mechanical equilibrium, or quasi- equilibrium. 

One may now ask whether the mechanical equilibrium is stable in the sense that it can 
last forever (assuming that the whole system remains perfectly isolated from the outside 
world) , or will the gases find ways to exchange energy through the piston and eventually 
bring the system to a thermal equilibrium? It was claimed in some textbooks, based 
on a simplistic interpretation of the laws of thermodynamics, see below, that indeed the 
mechanical equilibrium could persist "forever", cf. [G, GF] for some history. 

On the other hand. Landau and Lifshitz [LL], Feynman [F] and many others argued 
intuitively that the system should converge from the mechanical equilibrium to a thermal 
equilibrium. They predicted that the cooler compartment should gradually heat up and 
the hotter one cool down, while the piston slowly moves from the cooler side to the hotter 
side, so that the pressure balance is maintained until the temperatures are equalized, and 
the piston makes its final stop. 

The confusion about the evolution of the gas after the estabhshment of mechanical 
equilibrium is due to the following: Heat conduction through a wall is normally associated 
with the internal motion of the molecules of the wall colliding with those of the gas 
and thus exchanging momentum and energy. However, the piston and the walls in our 
idealized model are supposed to be rigid, solid, and structureless bodies and the gas 
atoms bounce off them elastically. This idealization is exactly the reason why the gases 
in the different compartments could be in equilibrium at different temperatures when 
the piston was clamped. The undamped piston, on the other hand, interacts with gas 
atoms as a whole, i.e. as one huge and massive molecule. It then makes tiny microscopic 
movements (vibrations) induced by collisions with atoms on both sides. Hence, some 
microscopic exchange of momentum and energy does take place. But these microscopic 
vibrations of the piston are not part of macroscopic thermodynamics, in which the action 
of the piston on the gas in each compartment is regarded as an external mechanical force. 
Under this condition (and assuming that the piston has no entropy of its own) the second 
law of thermodynamics would predict that the entropy of the gas, as it goes from some 
initial equilibrium state to a final equilibrium state, could not decrease. When the gas in 
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each compartment is in equilibrium, its thermodynamic entropy is known to be [LL, Ca] 

Si = Ni [- logP, + (1 + 3/2) logT,] + f{Ni), t = L,R 

where Nl = ulVl and Nji = ujiVr denote the number of atoms in the gases, and the 
exphcit form of f{Ni) is irrelevant for us, since its value docs not change in time. Now, 
if our system does evolve from mechanical equilibrium to thermal equilibrium, keeping 
the pressure balance Pl = Pr and the total kinetic energy ^kBiNiTi + N^Tji) fixed, 
then one can easily compute (we leave this as an exercise) that the pressure of the 
gases stays constant in time and the total entropy of the system S = Sl + Sr grows 
until it reaches its maximum at the point of thermal equilibrium. At the same time, 
the entropy Sr decreases, while Sl increases, since Tr goes down, Tl goes up, and the 
pressure Pr — Pl remains constant. This decrease of the entropy would, as already 
noted, violate the second law of thermodynamics, if the piston remained mechanical, see 
further discussions in [Li, CPS] and critical remarks in [G, GF]. 

Therefore, the evolution of the system beyond the mechanical equilibrium cannot be 
described by macroscopic thermodynamics (beyond the statement that any evolution in 
an isolated macroscopic system will not decrease the total entropy) . The actual evolution 
is a result of microscopic energy transfer between the gases via collisions with the piston. 
This process is purely microscopic and, in a sense, counterintuitive, as we explain next. 
Under the collisions with gas atoms on both sides the piston vibrates, i.e. it jiggles back 
and forth. When the piston moves toward the hotter side, the atoms of the hotter gas 
bounce off the piston with an increased speed and so gain energy, while the atoms of 
the cooler gas collide with the piston and slow down, hence lose some energy. When the 
piston moves toward the cooler side it is vice versa. Since, on the average, the hotter 
gas must cool down and the cooler gas must heat up, one may conclude that the piston's 
movements toward the cooler side dominate. On the other hand, the piston has to slowly 
move toward the hotter side in order to maintain the pressure balance, sec above, so its 
displacements in the direction of the hotter gas actually dominate. It is not quite clear 
how these seemingly opposite trends manage to coexist. Some physicists joke about a 
"conspiracy" between the microscopic vibrations of the piston and the incoming atoms 
of the gases [GF, GP]. In the words of Gallen [Ga], "the movable adiabatic wall presents 
a unique problem with subtleties" . 

In order to understand the mechanism of the heat transfer across the piston at me- 
chanical equilibrium [P^ = Pr), one usually considers the simplest gas of noninteracting 
particles, that is an ideal gas. As early as in 1959, Lebowitz [LI] studied a piston inter- 
acting with two infinite reservoirs filled with ideal gases held at different temperatures 
Tl 7^ Tr. His piston also interacted with an external potential, e.g. a spring, and it could 
therefore come to a stationary nonequilibrium state under the influence of the infinite 
reservoirs. He used an approximation by a Markov process and found the distribution 
of the piston velocity to be Maxwellian corresponding to some intermediate tempera- 
ture T G {Tl,Tr), which led to a systematic heat transfer between the gases. Recently, 
Gruber and others [GF, PG, GP] used kinetic theory to study a freely movable piston 
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of mass M interacting with two infinite ideal gases of atoms of mass m <S M at equal 
pressures but different temperatures. They use the expansion of the Boltzmann equation 
in e = m/M to show that a macroscopic heat flux across the piston does occur whenever 
Tl 7^ Tr, hence the system gradually approaches thermal equilibrium. They also found 
a stationary distribution of the piston velocity, whose average value is given by 



(it is independent of the gas densities). We note that if < Tr, then (V) > 0, confirming 
our previous observation that the piston moves from the cooler side to the hotter side. 

Equation (1.1) shows that the average velocity of the piston is different from zero, 
albeit just of order 0{m/M), despite the perfect pressure balance Pl = Pr. We note, 
however, that for a macroscopic-size piston the ratio m/M is so small that the time it 
takes the piston to cover any noticeable distance is much longer than the age of the 
universe [GF], so such a phenomenon cannot be observed experimentally. 

We conclude that the evolution of the system proceeds in two different stages. The 
first one is the convergence to a mechanical equilibrium, which is relatively fast and can, 
in principle, be computed on the basis of macroscopic equations. The second stage is 
the transition of the system from mechanical equilibrium to thermal equilibrium. This 
process is very slow and much less understood. 

In addition, for the ideal gas (in which the atoms do not interact) a new problem arises. 
At time t = 0, before the piston is released, the gas atoms move independently of each 
other - every atom bounces off the walls and the clamped piston, without exchanging 
momentum or energy with other atoms. Therefore, the velocity distribution does not 
have to be Maxwellian. A stationary state of the ideal gas can be described by any 
Poisson process with a uniform spacial density and a symmetric velocity distribution. 
For example, half of the atoms may move toward the piston with unit velocity v = 1 
and the other half - in the opposite direction with velocity v = —1, and this state will 
be stationary. However, once the piston is released, the atoms start interacting with 
each other, indirectly, via collisions with the piston. This provides a way to exchange 
momentum and energy between the atoms. One can expect that these interactions will 
lead, ultimately, to a true thermal equilibration, when the velocity distribution becomes 
Maxwellian, as we explain in Section 5. This process, however, may take even longer 
than the equilibration of the mean kinetic energies described above. 

To consider this new process in its "pure" form, we assume that initially the sys- 
tem is already in a homogeneous state - the gas density is constant across the entire 
cylinder and the velocity distribution is the same in both compartments (but different 
from Maxwellian). Then there seem to be no forces of any kind that would drive the 
piston anywhere. In particular, when the piston is initially placed in the middle of the 
cylinder, then by symmetry there should be no reason for it to move either way! On 
the other hand, the system is not in equilibrium until the velocity distribution becomes 
Maxwellian, hence it should find ways to evolve toward equilibrium, thus changing its 
macroscopic state. We discuss this further in Section 5. 




(1.1) 
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One can also consider a simpler case when that the container is infinitely long on both 
sides of the piston and the ideal gases have infinite number of atoms, as in [LI, GP, PG]. 
In that case the problem reduces to the classical Rayleigh gas - a big massive particle 
submerged in an ideal gas. In particular, our piston in an infinite cylinder becomes a 
one-dimensional Rayleigh gas, which we describe in some detail. 

Let a heavy tagged particle (called molecule) of mass M move on a fine under elastic 
collisions with atoms of mass m of an ideal gas with a uniform density n and some velocity 
distribution f{v) dv. Denote by X{t) and V{t) = X{t) the position and velocity of the 
molecule at time t. Even though f{v) need not be Maxwellian, the velocity function V{t) 
and the coordinate function X{t) can be approximated by certain Gaussian stochastic 
processes: 

Theorem 1.1 (Holley [H]) Let the density f{v) he symmetric f{v) = f{~v) and have 
a finite fourth moment J v'^fiy) dv < oo. Then for every finite to < oo, the function 
V{t)y/M on the interval [0,to] converges, in distribution, as M,n — > oo and M/n — > 
const, to an Ornstein-Uhlenbeck velocity process Vt, while X(t)y/M converges to an 
Ornstein-Uhlenbeck position process Xt- 

An Ornstein-Uhlenbeck process {X^, Vt) is defined by [Ne] 

dXt = Vt dt, dVt = -aVt dt + \fD dWt 

where a > 0, > are constants and Wt a Wiener process. The Ornstein-Uhlenbeck 
position process Xt converges in an appropriate limit (e.g. a — > oo, a^/D = const) to a 
Wiener process. 

We note that the typical velocity of the molecule V{t) is of order 0(1/\/M), which 
agrees with the equipartition of energy in the system requiring that average energies of 
all particles be equal, i.e. M{V'^) = m{v^) = m J v^f{v) dv. 

Diirr et al [DGL] extended the above theorem to arbitrary dimension and to asymmet- 
ric velocity distributions. The main technical difficulty in the proof of this theorem comes 
from the so called recollisions, which occur when an atom collides with the molecule more 
than once. Recollisions result in intricate autocorrelations in the process {X{t),V{t)), 
which otherwise would be Markovian. The proof essentially consists in estimating the 
undesirable effect of recollisions and showing that it vanishes in the limit M — > oo. 

When the gas is confined in a finite cylinder, though, the effect of recollisions becomes 
crucial. All atoms will travel to the walls, bounce off it and come back to the piston for 
more and more collisions. The induced autocorrelations will build up. There is no 
standard techniques available to estimate (let alone eliminate) the effect of recollisions 
in general, but we make partial progress in this direction, see Section 4. 

In summary, the piston problem raises serious mathematical questions and even leads 
to confusions in the physical theories. The "notorious piston", as it is known among 
physicists, again attracted much attention recently due to a series of papers [GF, GP, 
PG, LPS] where a more extensive mathematical apparatus was developed. At the same 
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time, many new numerical experiments led to better theoretical understanding of the 
underlying dynamics but also raised some new questions. We emphasize that very few 
rigorous results are available, even for ideal gases, apart from the Rayleigh-type stochastic 
approximations in the infinite cylinder mentioned above. 

We study the piston in a finite cyhnder filled with ideal gases. Since our gases are 
ideal, we will not need to assume that the velocity distribution of atoms is Maxwellian. 
Since autocorrelations induced by recoUisions present a major difficulty, we specify the 
initial state in such a way, that during a certain interval of time each gas atom collides 
with the piston at most twice. The main goal of our work is to describe rigorously 
the dynamics of the piston during that time interval. We show that, in an appropriate 
limit, the evolution of the piston and the gas converges to a deterministic process, which 
satisfies a certain closed system of differential equations. The assumptions that we make 
here simplify technical considerations but by no means reduce the problem to a triviality. 
In fact, many intriguing questions still remain open in our context, and we discuss them 
in the last two sections of the paper. 

Precise statement of problem and main results. Consider a cubical domain of 
size L separated into two parts by a movable wall (piston). Each part of A/, contains a 
gas of noninteracting particles of mass m—1. The particles collide with the outer (fixed) 
walls of Ax, and with the moving piston elastically. The piston has mass M = and 
moves along the x-axis under the collisions with the gas particles on both sides. The size 
L of the cube is a large parameter of our model, and we are interested in the behavior as 
L — > oo. We will assume that is proportional to the area of the piston, i.e. Ml ~ -L^, 
and the number of gas particles N is proportional to the volume of the cube A^,, i.e. 
N L^, while the particle velocities remain of order one. 

The position of the piston at time t is specified by a single coordinate X = Xx(t), 
< X < L, its velocity is then given by ^ = Vi,(t) = X^it). Since the components of 
the particle velocities perpendicular to the x-axis play no role in the dynamics, we may 
assume that each particle has only one coordinate, x, and one component of velocity, v, 
directed along the x-axis. 

When a particle with velocity v collides with the piston with velocity V, their velocities 
after the collision, v' and V, respectively, are given by 

V'^(l-s)V + ev (1.2) 

v' = -{1 - e)v + {2 ~ e)V (1.3) 

where e = 2m/ {M + m). We assume that M + m = 2mL^/a, where a > is a constant, 
so that 

_ 2m _ 

^'M + m^L^ ^ ' 

When a particle collides with a wall at a; = or a; = L, its velocity just changes sign. 

The evolution of the system is then completely deterministic, but one needs to specify 
the initial conditions. We shall assume that the piston starts at the midpoint Xl{0) = 
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L/2 with zero velocity Vl(0) = (see also Section 2). The initial configuration of gas 
particles and their velocities is chosen at random as a realization of a (two-dimensional) 
Poisson process on the (a;,f)-plane (restricted to < a; < L) with density L'^pl{x,v), 
where Pl{x, v) is a function satisfying certain conditions, see below, and the factor of is 
the cross-sectional area of the container. This means that for any domain D C [0, L] x 
the number Nd of gas particles (x, v) E D at time t — has a Poisson distribution with 
parameter 

Xd — jj Pl{x, v) dxdv 

For any two nonoverlapping domains, say Di fl D2 = 0, the corresponding numbers N^j^ 
and are statistically independent. We remark that the total number of gas particles 
is a Poisson random variable, too. The total energy and the total initial momentum 

are random as well. 

Let VLl denote the space of all possible configurations of gas particles in (i.e., all 
countable subsets of [0, L] x R^). For each reahzation^ a; e the deterministic piston 
trajectory will be denoted by XL{t,ijj) and its velocity by VL{t,uj). 

The above model is a mechanical system whose dynamical characteristics Xiit.uj) 
and Viit^u)) depend on the large parameter L and, for each L, are random (depend on 
a;). 

In order to obtain a deterministic description of the dynamics of the piston one needs 
to take a limit as L — > 00 and simultaneously rescale space and time. We introduce new 
space and time coordinates by 

y = x/L and T — t/L. (1.5) 

which corresponds to Euler scaling for the hydrodynamical limit transition. We call y 
and r the macroscopic ( "slow" ) variables, as opposed to the original microscopic ( "fast" ) 
X and t. Now let 

Yl{t,uj) = XL{rL,uj)/L, Wl{t,u;) = Vl{tL,uj) (1.6) 

denote the position and velocity of the piston in the macroscopic context. The initial 
conditions are then ^^(0) = Xi(0)/L = 0.5 and 1^(0) = V{0) = 0. 

It is now very natural to assume that the initial density pl{x,v) agrees with our 
rescaling: 

PLix, v) = 7ro(x/L, v) (1.7) 

where the function 7ro(|/, v) is independent of L. Without loss of generality, we can assume 
that TTo is normalized so that 

r-l roo 



/ / ■Ko{y,v)dvdy 

Jo J-00 



^Technically, it is possible that two or more particles collide with the piston simultaneously, and then 
the dynamics will no longer be defined, but multiple collisions are known to occur with probability zero 
[H], so we will ignore such anomalies. 
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Then the mean number of particles in the entire container A^, is exactly equal to L^: 

E{N) ^ jj L'^pl{x, v) dv dx = 

where E{-) is the expected value. 

Furthermore, we assume that the function 7ro{y,v) satisfies several technical require- 
ments stated below. The meaning and purpose of these assumptions will become clear 
later. 

(PI) Smoothness. 7io{y,v) is a piecewise function with uniformly bounded partial 
derivatives, i.e. Idno/dy] < Di and \dT^Q/dv\ < Di for some Di > 0. 

(P2) Discontinuity lines. T:o{y,v) may be discontinuous on the hne y = Yl{0) (i.e., "on 
the piston"). In addition, it may have a finite number (< Ki) of other discontinuity 
lines in the {y, i')-plane with strictly positive slopes (each line is given by an equation 
V — f{y) where f{y) is and < ci < f'{y) < C2 < oo). 

(P3) Density bounds. Let 

no{y,v) > TTmin > for Vi < \v\ < V2 (1.8) 

for some < Vi < V2 < 00, and 

SUpTlo{y,v) = TTmax < OO (1.9) 

The requirements (1.8) and (1.9) basically mean that 7io{y,v) takes values of order 
one. 

(P4) Velocity "cutoff". Let 

7ro(?/,f)=0, if \v\<v^in or \v\>Vma^ (1-10) 

with some < Vmin < ^max < OO. This means that the speed of gas particles is 
bounded from above by fmax and from below by fmin- 

(P5) Approximate pressure balance. 7ro{y, v) must be nearly symmetric about the piston, 
i.e. 

\Tro{y,v)-no{l-y,-v)\<£o (1.11) 
for all < 1/ < 1 and some sufficiently small £o > 0. 

The requirements (P4) and (P5) are crucial. We will see that they are made to ensure 
that the speed of the piston \ VL{t, uj)\ will be smaller than the minimum speed of the gas 
particles, with probability close to one, for times t = 0{L). Such assumptions were first 
made in [LPS]. 
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We think of A, i^i, (P1)-(P4) as fixed 

(global) constants and in (P5) as an adjustable small parameter. We will assume 
throughout the paper that Sq is small enough, meaning that 

£0 < £o(-Dl, Ki, Cl, C2, t"!, f2, t-min, ^'max, TTmin, TTmax) 

It is important to note that the hydrodynamic limit does not require that Eq — >■ 0. The 
parameter £0 stays positive and fixed as L — > (X). 
Now we state our main result: 

Theorem 1.2 There is an L-independent function Y{t) defined for all r > and a 
positive r* ~ 2/ v^ax ( actually, r* — >• 2/ v^ax as Eq ^ ), such that 

sup \Yl{t,uj)-Y{t)\^Q (1.12) 

0<r<r, 

and 

sup \Wl{t,uj)-W{t)\^Q (1.13) 

0<T<T, 

in probability, as L ^ 00. Here W{t) — Y{t). 

This theorem establishes the convergence in probability of the random functions 
Yl{t,uj)^ W{t^uj) characterizing the mechanical evolution of the piston to the determin- 
istic functions Y{t)^ W{t), in the hydrodynamical limit L — > cxd. 

The functions Y{t) and W{t) satisfy certain (Euler-type) differential equations stated 
in the next section. Those equations have solutions for all r > 0, but we can only 
guarantee the convergence (1.12) and (1.13) for r < r*. What happens for r > r*, 
especially as r — > cxo, remains an open problem. Some numerical results and heuristic 
observations in this direction are presented in Section 5. 

Remarks. The function Y{t) is at least and, furthermore, picccwise . On the 
interval (0,r*), its first derivative W = Y (velocity) and its second derivative A = Y 
(acceleration) remain £o-sniall: sup^|iy(T)| <const-£o and sup^|yl(r)| <const-£o, see 
the next section. 

We will also estimate the speed of convergence in (1.12) and (1.13). Precisely, we 
show that there is a tq > (tq ~ 1/t'max) such that 

\YL{T,u;)-Y{T)\ = 0{lnL/L) 

for < T < To and 

\YL{T,co)-Yir)\=Oi\nL/L'/') 

for To < T < r*. The same bounds are valid for \Wl{t,uj) — W{r)\, see Sections 3 and 
4. These estimates hold with "overwhelming" probability, specifically they hold for all 
uenicflL such that P{ni) = 1 - 0(L-''^^). 
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2 Hydro dynamical equations 

The equations describing the deterministic function Y[r) involve another deterministic 
function - the scaled density of the gas TT{y, v, r). Initially, 7r(y, v, 0) = 7ro(y, v), and for 
T > the density Tr{y, v, r) evolves according to the following rules. 

(HI) Free motion. Inside the container the density satisfies the standard continuity 
equation for a noninteracting particle system without external forces: 

for all y except y — 0, y — 1 and y — Y{t). 

Equation (2.1) has a simple solution 

7r{y,v,T) ^ 7r{y -vs,v,T - s) (2.2) 

for < s < r such that y — vr ^ {0, Y{t — r), 1} for all r G (0, s). Equation (2.2) has one 
advantage over (2.1): it applies to all points {y,v), including those where the function tt 
is not differentiable. 

(H2) Collisions with the walls. At the walls y — and y — 1 we have 

tt{0,v,t) =7r{0,-v,T) (2.3) 
7r(l,^;,T) = 7r(l,-^;,T) (2.4) 

(H3) Collisions with the piston. At the piston y — Y{t) we have 

7i{Y{t) -0,v,t) = 7fiY{T) -0,2W{t) -v,t) for v<W{t) 
tt{Y{t) + Q,v,t) = ti{Y{t) + Q,2W{t) -v,t) for v > VF(t) (2.5) 

where v represents the velocity after the collision and 2W{t) — v that before the 
collision; here 

= ^Y{r) (2.6) 
is the (deterministic) velocity of the piston. 

It remains to describe the evolution of W{t). Suppose the piston's position at time 
r is y and its velocity W. The piston is affected by the particles {y, v) hitting it from 
the right (such that y = Y + and v < W) and from the left (such that y = Y — and 
V > W). 
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(H4) Piston's velocity. The velocity W — W{r) of the piston must satisfy the equation 
/ {v-Wf7r{Y -0,v,T)dv^ / {v - WfTr{Y + 0,v,t) dv (2.7) 

see also an additional requirement (H4') below. 

In physical terms, (2.7) is a pressure balance: the piston "chooses" velocity W so that 
the pressure of the incoming particles balances out. Equation (2.7) is instrumental for 
our deterministic approximation of the piston dynamics. 

One can combine the two integrals in (2.7) into one by introducing the density of the 
particles colliding with the piston ( "density on the piston" ) by 

, , f 7r(F + 0,w,r) ii v<W 

^("'^;^'^) = |^(F-0,.,r) v>W (2.8) 

Then (2.7) can be rewritten as 

' {v - W{T) f sgn(y - W{t)) q{v, r; F(r), W{t)) dv = 



We also remark that for r > 0, when (2.5) holds, 

/ V7i(Y — 0, v, r) dv J V7i(Y + 0, v, r) dv 



W(t) = 



J 7i(Y - 0,v,t) dv JniY + 0,v,T)dv 



i.e. the piston's velocity is the average of the nearby particle velocities on each side. 

The system of (hydrodynamical) equations given in (H1)-(H4) is closed and, given 
appropriate initial conditions, should completely determine the functions Y{t), W{t) 
and 7r{y, v, r) for r > 0, as we will see shortly. 

To specify the initial conditions, we set 7r{y,v,0) = 7ro{y,v) and Y{0) = 0.5. The 
initial velocity W{0) does not have to be specified, it comes "for free" as the solution of 
the equation (2.7) at time r = 0. It is easy to check that the initial speed |Vr(0)| will be 
smaller than Vmin, in fact W{0) — > as £o — in (P5). 

We first determine conditions under which equation (2.7) has a solution W. Let 

^sup(^) = sup{w : 7r{Y - 0, v, r) > 0} 
(with the convention that the supremum of an empty set is — oo) and 

v+f (r) = mi{v : 7r{Y + 0, v, r) > 0} 
(similarly, the infimum of an empty set must be set to +oo). 
Lemma 2.1 We have three cases: 
(a) // ^sup > t'inf or Wg^p = fjnf e IR, then (2.7) has a unique solution W G [fjnf, I'sup]. 
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(b) IfVg^p < v^f, then the solutions of (2.7) occupy the entire interval ['t'sup; '^infl- 

(c) If v~^p = fj^f = oo or = = —oo, then (2.7) has no real solutions. 

Proof. In the case (a), the difference between the left hand side and the right hand side of 
(2.7) is a continuous and strictly monotonically decreasing function of W. For W < 
it is positive, and for W > v'^^ negative. The rest of the proof goes by direct inspection. 
□ 

It is easy to show (we do not elaborate) that under our assumptions (P1)-(P4) for 
every r > the density 7r{y, v, r) has a compact support on the y, v plane, i.e. 7r{y, v, r) = 
for all \v\ > fmax(''")- Therefore, the "no solution" case (c) never occurs. The multiple 
solution case (b) is very unlikely, but not impossible. If that happens, the velocity W{t) 
must be defined uniquely by an additional requirement: 

(H4') If W{t - 0) e fcp,^'^], we define W{t) by continuity, W{t) = W{t - 0). If 
W{t — 0) < v~p or W{t — 0) > v^f, we set W{t) = v~^p ov W = v^^, respectively. 

This completes the definition of W{t) started by (H4). 

For generic piecewise smooth densities Tr{y, v, r), the velocity W{t) is continuous, but 
in some cases the continuity of W{t) might be broken. The following simple lemma will 
be helpful, though: 

Lemma 2.2 Suppose that for every r e [a, b] the density 7r{y, v, r) is piecewise and 
has a finite number ofC^ smooth discontinuity lines on the y, v plane with positive slopes, 
as we require of -KQ{y,v) in Section 1. Then W{t) will be continuous and piecewise 
differentiable on the interval [a, b] . 

We now pause to make a few remarks. The piston mass is never used in our equations, 
because its macroscopic mass is zero. Indeed, for the mechanical system described in 
Section 1, the piston mass is ~ L^, while the total mass of the gas particles is ~ L^, 
hence the relative mass of the piston vanishes as L — > oo. Consider now the total 
(macroscopic) mass of the gas 




and the mass in the left and right compartments, separately. 





and the total kinetic energy 




The following lemma is left as a (simple) exercise: 
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Lemma 2.3 The quantities Attot? M.L, M.R, and £tot remain constant in r. 



The main equation (2.7) also preserves the total momentum / / v7r{y,v,T) dvdy, but 
this quantity changes due to collisions with the walls. 

Remark. Previously, Lebowitz, Piasecki and Sinai [LPS] studied the piston dynamics 
under essentially the same initial conditions as our (P1)-(P5). They argued heuristically 
that the piston dynamics could be approximated by certain deterministic equations in the 
original (microscopic) variables x and t. In fact, the present work grew as a continuation 
of [LPS]. The deterministic equations found in [LPS] correspond to our (2.2)-(2.6) with 
obvious transformation back to the variables x,t, but our main equation (2.7) has a 
different counterpart in the context of [LPS], which reads 



/ {v- V{t)f7r{Y - 0, V, t)dv- {v - V{t)fn{Y + 0, v, t) dv 

Jv J —oo 



(2.9) 



Here X — X{t) and V — V{t) — X{t) denote the deterministic position and velocity of 
the piston and T^{x,v,t) the density of the gas (the constant a appeared in (1.4)). We 
refer to [LPS] for more details and a heuristic derivation of (2.9). Since (2.9), unlike our 
(2.7), is a differential equation, the initial velocity V{Q) has to be specified separately, 
and it is customary to set V"(0) = 0. Equation (2.9) can be reduced to (2.7) in the limit 
L — > oo as follows. One can show (we omit details) that (2.9) is a dissipative equation 
whose solution with any (small enough) initial condition ^(0) converges to the solution 
of (2.7) during a t-time interval of length ~ InL. That interval has length ~ L~^lnL 
on the T axis, and so it vanishes as L oo, this is why we replace (2.9) with (2.7) and 
ignore the initial condition V(0) when working with the thermodynamic variables r and 
y. For the same reasons, it will be convenient to reset the initial value of the piston 
velocity in the mechanical model of Section 1 to from y(0) = to V{Qi) — VF(0), see 
Theorem 3.5 below. The equation (2.9) will not be used anymore in this paper. 

We now describe the solution of the hydrodynamical equations (H1)-(H4) in more 
detail. Assume that for some r > the gas density 7r{y,v,T) satisfies the following 
requirements, similar to (P1)-(P5) imposed on the initial function 7ro(y,t') in Section 1: 

(PI') Smoothness. 7r{y,v,T) is a piecewise function with uniformly bounded partial 
derivatives, i.e. \d7r/dy\ < D[ and \d7r/dv\ < D[ for some D[ > 0. 

(P2') Discontinuity lines. 7T{y,v,T) has a finite number (< K[) of discontinuity lines 
in the {y, i')-plane with strictly positive slopes (each line is given by an equation 
V = f{y) where f{y) is and < c'l < f'{y) < 4 < oo). 

(P3') Density bounds. Let 

7r(t/,'t;,r) > > for < [-ul < t;^ (2.10) 
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for some < v[ < V2 < 00, and 



sup7r(7/,i;,T) = 7r;„3^ < 00 (2.11) 



(P4') Velocity "cutoff". Let 

7r(y,v,r)^0, if \v\ < v'^^ or |^;| > (2.12) 

with some < v'^^ < v'^^ < 00. 

Lastly, we want to assume, similarly to (P5), that Tr{y,v,T) is nearly symmetric about 
the piston, but this assumption requires a little extra work, since the piston docs not 
have to stay at the middle point Y{0) = 0.5. For every Y G (0, 1) denote by hy the 
unique homeomorphism of [0, 1] such that /iy(0) = 1, /ty(l) = 0, /iy(F) = Y and hy is 
linear on the subsegments [0, F] and [Y, 1]. Next, we consider [0,F] x R as a manifold 
in which points (0, v) and (0, —v) are identified for all v > 0, and so are the points {Y, v) 
and (y, —v) for v > 0. Similarly, let [Y, 1] x IR be a manifold in which one identifies (1, v) 
with (1, —v) and (1^,^) and (Y, —v) for all v > 0. We denote by dy the distance on each 
of these two manifolds induced by the Euclidean metric {dy^ + dv'^Y^'^. The reason why 
we need this special distance will be clear later, in the proof of Proposition 2.10. 

(P5') Approximate pressure balance. We require that 

|r(r) -0.5| < e'o (2.13) 

and for any point {y, v) with < y < 1 and < |^| < ^max there is another point 
(y*,f*) "across the piston", i.e. such that {y — Y){y^ — F) < 0, where Y = F(r), 
satisfying 

dy{{y.,v.),{hY{y),-v))<e'^ (2.14) 

and 

k(y,^;,T) - 7r(y*,^;*,r)| < £0 (2.15) 
for some sufficiently small £q > 0. In addition, we require that 



4 < C^£o (2.16) 



with some constant Cq > 0. 



Actually, the map (t/,f) 1— > {y*,v*) involved in (P5'), which we will denote by Rr, is 
one-to-one and will be explicitly constructed below, in the proof of Proposition 2.10. 

Again, we think of D[, K[, c[, 4, v[, v'^, v'^,^, v'^^^, 'n'^,^, and now also Cq, as 

global constants. They must be bounded on the time interval on which we consider the 
dynamics (and I'^in? "^min must be bounded away from zero), hence we may treat all these 
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constants as independent of r. By (2.16), e'q is, just like Eq in (P5), a small adjustable 
parameter. 

Now we derive rather elementary but important consequences of the above assump- 
tions. Since the density 7T{y, v, r) vanishes for \v\ < v'^^^, so does the function q{v, r; Y, W) 
defined by (2.8). Moreover, for all \W\ < v'^^^, the function q{v, r; Y, W) will be indepen- 
dent of W and can be redefined by 

g(.,r;y) = | + " < ° (2.17) 

Also, the equation (2.7) can be simplified: the factor sgn(t' — can be replaced by 
sgnw. Then, expanding the squares in (2.7) reduces it to a simple quadratic equation for 
W: 

QoW^' - 2Qiiy + Q2 = (2.18) 

where 

Qo = / %Z'^v-q{v,T;Y)dv (2.19) 

Qi^ j vsgnvq{v,T;Y)dv (2.20) 

Q2 = j v^sgnvq{v,T;Y)dv (2.21) 
with Y — Y{t). The integrals Qq, Qi, Q2 have the following physical meaning: 

mQo = mL-mR 

mQi ^Pl-Pr 
mQ2 = 2{eL - br) 

where mL-iPLie-L represent the total mass, momentum and energy of the incoming gas 
particles (per unit length) on the left hand side of the piston, and mR,pR,eR - those 
on the right hand side of it. The value Q2 also represents the net pressure exerted on 
the piston by the gas if the piston did not move. Of course, if Q2{t) = 0, then we 
must have W{t) — 0, which agrees with (2.18). The following lemma easily follows from 
(P1')-(P5'). It means that the function q{v,T]Y{T)) is nearly symmetric in v about 
^; = 0. 



Lemma 2.4 For any smooth function f{v) defined for v > we have 



f{v)q{v,T;Y{T))dv - / f{-v)q{v,T;Y{T))dv 
J-00 

where the factor C/ > depends on f but not on Eq- 



<CfEo 
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Convention. We call constants that do not depend on our small adjustable parameter 

£q involved in (P5) and (P5') global constants (such as Cf in the above lemma). All the 
constants in the requirements (P1)-(P5) and (Pr)-(P5') are global, except Sq itself and 
the related Eq. In many cases, we will denote various global constants by Q, i > 0, or 
just by C. 

Lemma 2.4 implies that Qo and Q2 are small, more precisely 

max{|Qo|,|Q2|} <C£o (2.22) 

where C > is a global constant. At the same time, the assumption (P3') guarantees 
that 

Ql > Ql,min > (2.23) 

where Qi,mm is another global constant. 

If Eq is small enough, there is a unique root of the quadratic polynomial (2.18) on the 
interval (— "i^min) "^min); which corresponds to the only solution of (2.7). Since this root is 
smaller, in absolute value, than the other root of (2.18), it can expressed by 

^^e.-yep^a (2.24) 

Qo 

where the sign before the radical is "— ", not "+". Of course, (2.24) apphes whenever 
Qo 7^ 0, while for Qo = we simply have 

= ^ (2.26) 

Corollary 2.5 If Eq is small enough, then 

\W{T)\<Beo<v'^j3 (2.26) 

with some global constant B > 0. 

Proof. This immediately follows from equations (2.22)-(2.25). □ 

We now make an important remark. 

Remark (Extension). Consider the density of the incoming gas particles on the left 
hand side of the piston, i.e. 7r{y,v,T) for y = Y{t) — and v > v^nm- This function 
"terminates" on the piston, i.e. has a discontinuity in y at y = Y{t). But it can be 
naturally extended smoothly "across the piston", i.e. for y > Y{t) if one ignores the 
interaction of the gas coming from the left compartment with the piston at times s e 
(r — S,t) and applies the rule (HI) instead, as if the gas "passed through the piston". 
This defines a smooth extension of 7r(?/,f,r) from the region y < Y{t) to the region 
Y{t) < y < Y{t) + 0{S) for all v > f'^in- This extension allows us to differentiate 
q{v,T;Y) defined by (2.17) with respect to Y for any v > v'^^^. A similar extension can 
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be made for the density Tr{y, v, r) from the region y > Y{t) to the region Y{t) — 0{S) < 
y < Y{t) for all v < — fmin, hence q{v,T;Y) becomes differentiable with respect to Y for 

V < ^"^min- We note that our extension can be unambiguously defined because we only 
need it for \v\ > f^j^ while the piston's velocity remains smaller than fniin- 

Now the quantities Qq, Qi, and Q2 defined by (2.19)-(2.21) become differentiable in 

Y for each fixed r, and the assumptions (P1')-(P4') easily imply that 

\dQi/dY\<Ci, i = 0,l,2 (2.27) 
where Ci > is a global constant. 

Corolleiry 2.6 The piston acceleration A{t) = dW{T)/dT satisfies 

\A{t)\ < Ceo (2.28) 

with a global constant C > 0. 

Proof. We differentiate the quadratic equation (2.18) with respect to r and get 

. ^ {dQo/dT)W^ - 2{dQi/dr)W + {dQ2/dT) 
~ 2(Qi ~ Q^W) 

Clearly, the denominator is bounded away from zero, and the numerator has an upper 
bound of order sq, because \dQi/dT\ = \{dQi/dY)W\ < const-^o by (2.27) and (2.26). □ 

More importantly, we can now derive the existence and uniqueness of the solution of 
the hydrodynamical equations (H1)-(H4) as long as the conditions (Pr)-(P5') continue 
holding: 

Lemma 2.7 // the hydrodynamical equations (H1)-(H4) have a solution on an interval 
< T < T and the conditions (P1')-(P5') hold on this interval, then the solution is 
unique at t — T and can be extended immediately beyond the point t — T. 

Proof. The only differential equation in our system (H1)-(H4) is (2.6), in which W{t) 
is the root of the quadratic equation (2.18) given by (2.24)-(2.25). Due to the above 
Exteusiou Remark we can think of W as an implicit function of Y, i.e. effectively W — 
F{Y,t). Then the differential equation (2.6) takes a canonical form 

^F(r) = F(F(r),r) (2.29) 

For this equation to have a Tinique solution, it suffices that F{Y, r) has a bounded partial 
derivative with respect to Y. 

Since W is a root of the quadratic equation (2.18), we can differentiate (2.18) with 
respect to Y and get 

dF{Y,T) _ {dQo/dY)W^ - 2{dQi/dY)W + {dQ2/dY) 
dY ~ 2(Qi - QoW) 
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We already know that the denominator is bounded away from zero. It follows from (2.27) 
that the numerator stays bounded above, hence 

dF{Y, t) 



dY 



< K (2.30) 



with a global constant k > 0. □ 

Next we consider the evolution of a point {y, v) in the domain 

g:={{y,v): < y < 1} 

under the rules (H1)-(H3), i.e. as it moves freely with constant velocity and collides 
elastically with the walls and the piston. Denote by {y-,-, Vr) its position and velocity at 
time r > 0. Then (HI) translates into ijr = Vr and v-,- — whenever y^ ^ {0, 1,Y{t)}, 
(H2) becomes {yr+o,Vr+o) — {Ut-o, —Vt-q) whenever y^_o G {0, 1}, and (H3) gives 

(yr+o, Vr+o) = (yr-0, 2W^(t) - Vr-o) (2.31) 

whenever y-r-o — Y{'t)- Note that (2.31) corresponds to a special case of the mechanical 
collision rules (1.2)-(1.3) with e = (equivalently, m — 0). Hence the point {y,v) moves 
in Q as if it was a gas particle with zero mass. 

The motion of points in [y, v) is described by a one-parameter family of transfor- 
mations J-"^ '■ Q ^ Q defined by jF^d/ci'o) = {yr-iVr) for r > 0. We will also write 
^'^{UtjVt) — {yo,vo)- According to (H1)-(H3), the density 7r{y,v,T) satisfies a simple 
equation 

7r(yr, Vr, r) = ■K{J^~''{yr, v^), 0) = 7ro(|/o, vo) (2.32) 

for all r > 0. Also, it is easy to see that for each r > the map J-'^ is one-to-one and 

preserves area, i.e. det \DJ-"^{y,v)\ = 1. 

Now, because of (P4), the initial density 7io{y,v) can only be positive in the region 

:= {{y,v) : < y < 1, t-min < |?;| < t;max} 

hence we will restrict ourselves to points {y,v) G only. At any time r > 0, the 
images of those points will be confined to the region Q~^{t) :— J-'^{Q^). In particular, 
7^(?/,^^,r) = for {y,v) ^ G+{t). 

We now make an important observation. If a point {yT,Vr) collides with a piston 
whose velocity is slow, |W^(t)| <^ \vr\, they cannot recollide too soon: the point must 
travel to a wall, bounce off it, and then travel back to the piston before it hits it again. 
This is quantified in the following lemma: 

Lemma 2.8 Let a point (y^jV^-) G G~^{t) collide with the piston, i.e. = Y{t). Then 
during the interval (r, r + A) with 

^_ l-2Bsor 

it cannot recollide with the piston, i.e. y^ ^ Y{s) for s G {t,t + A), provided (P1')~(P5') 
continue holding during this interval. 
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Proof. The point's speed after the coUision is at least v'^^^ — 2BeQ and at most v'^^^-\-2BeQ. 
The piston cannot "catch up" with it, since |VF(r)| < ^j'^j^ — 2BeQ by (2.26). So, the 
point travels to the wall, bounces off it, and travels back to the piston, and all that will 
take time 

A > 2D/«,, + We,) 

where D = min{F(r), 1 - F(r)} > 0.5 - BeqT. □ 

Therefore, as long as (P1')-(P4') hold, the collisions of each moving point (t/^,t>^) e 
^^(r) with the piston occur at well separated time moments, which allows us to effectively 
count them. For (x, v) e ^+ 

N{y, V, r) = #{s e (0, r) : = Y{s), v, 7^ W{s)} 

is the number of collisions of the point {y,v) with the piston during the interval (0, r). 
For each r > 0, we partition the region G^{t) into subregions 

Q^iT):={J'^{y,v): {y,v)eg+ & N{y,v,T) = n} 

so Qni^) is occupied by the points that at time r have experienced exactly n collisions 
with the piston during the interval (0,t). 

Now, for each n > 1 we define t„ > to be the first time when a point {y^, Vr) E G^{t) 
experiences its {n + l)-st collision with the piston, i.e. 

T„ = sup{r>0: ^;++i(t) = 0} 

In particular, ti > is the earliest time when a point {yT,Vr) G G~^{t) experiences its 
first recollision with the piston. Hence, no recollisions occur on the interval [0,ri), and 
we call it the zero-recollision interval. Similarly, on the interval {ti,T2) no more than 
one recollision with the piston is possible for any point, and we call it the one-recollision 
interval. 

The time moment mentioned in Theorem 1.2 is the earliest time when a point 
(^y^^Vr) G g^if) either experiences its third collision with the piston or has its second 
collision with the piston given that the first one occurred after ri. Hence, r* < T2, and 
actually is very close to T2, see the next lemma. 

Lemma 2.9 Let (Pr)-(P5') hold on the interval {0,n/vmax + for some n > 1 and 
S > 0. Then, for all sufficiently small Eq 

for all 1 < k < n, where C > is a global constant that may depend on n. Also, 

- 2/finax| < CEq 
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Proof. The necessary lower bounds on follow from Lemma 2.8. The necessary upper 
bounds are just as easy to obtain, we omit details. □ 

It is clear at this point that the hydrodynamical equations (H1)-(H4) will have a 
unique and "well behaved" solution as long as the conditions (P1')-(P5') continue holding 
with some small Eq. Our next goal is to show that this is indeed the case. 

Proposition 2.10 Let T > 0. // the initial density 7iQ{y,v) satisfies (P1)-(P5) and 
in (P5) is small enough (for the given T), then the conditions (P1')-(P5') will hold on 
the interval < r < T . 

Note: The corresponding global constants in (P1')-(P5') will depend on T as specified 
below. 

Proof. The main idea is to show that the restrictions on 7r(y, v, r) imposed by (P1)-(P5) 
at r = "deteriorate" very slowly, as time goes on, so that (Pr)-(P5') will continue 
holding ( "propagate" ) with the respective global constants slowly changing in time. 

We first note that as long as (P1')-(P5') hold, the number of collisions grows at most 
linearly in r, i.e. on any interval (0, r) on which (Pr)-(P5') hold, every moving point 
{ys,Vs) G S~^{s) experiences at most tv'^^^ + 1 collisions with the piston, if Eq is small 
enough, see Lemmas 2.8-2.9. Next, we examine the conditions (Pr)-(P5') individually 
and show that each of them should hold up to time T, provided that the others do. 

We start with (PI). Due to (HI) we have 

d7r{y,v,T + s) d7r{y — sv,v,t) 
dy dy 

and 

d7r{y,v,T + s) d7r{y — sv,v,t) d7r{y — sv,v,t) 

s ~ 



dv dv dy 

for all s > such that the moving point located at {y, v) at time t + s did not experience 
collisions with the piston during the interval (r, t + s). Thus, between collisions with the 
walls and the piston, the partial derivatives of TT[y,v,r) can grow at most linearly with 
r. Collisions with the walls could only change the sign of the derivatives of p, but not 
their absolute values. 

Now consider the effect of interactions with the piston. We evaluate the partial 
derivatives of n{y,v,T) at a point {y,v) after a collision with the piston at some earlier 
time s G (0,t). For simplicity, assume that there are no other collisions of the moving 
point {y, v) with the piston or the walls on the interval (s, r). Then s satisfies the equation 

Y(s) ^y-(T-s)v (2.33) 

Due to (H3) and (HI) we have 

TT{y,v,T) = n{y — {t — s)v,v, s + 0) 

= 7i{y - (r - s)v, 2W - v,s - 0) 

= 7r{y - (r - s)v - {s - so)i2W - v), 2W - v, so) 
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where so < s is any earlier time (that we consider fixed) and W — W{s) is the piston 
velocity at the time of collision. Let fjo — y — {t — s)v — {s — Sq){2W — v) and Vq — 2W — v. 
Then 



d'K{y,v,T) dTi{yo,Vo,SQ) 



dy dy 



ds ^, ,dW ,ds 
l + v—-2{s- So)— - {2W - v) — 
dy dy dy 



^ d7r{yo,vo,so) ^ dW 
dv dy 

Difi^erentiating (2.33) with respect to y gives 

dY ds ^ ds 
ds dy dy 

hence 

ds 1 



dy W 

Also, 



dW dW ds A 



dy ds dy W — v 
where A — A(s) is the piston acceleration at the time of collision. Now, as long as 
(P1')-(P5') hold, we have W = O{eo) and A = O{eo), hence ds/dy = -y-^ + O{eo) and 
so 

d7r{y,v,T) ^ d7r{yo,vo,so) ^ ^ 
dy dy 

In other words, the piston (due to its low speed and acceleration) acts almost as a wall, 
which only changes the sign of dn/dy. A similar calculation (we omit it) holds for the 
partial derivative with respect to v. 

Thus, as long as (Pr)-(P5') hold, the density n{y,v,T) remains piecewise and its 
partial derivatives can grow at most linearly with r. 

Next, we check the condition (P2'). We begin with three special discontinuity lines 
that do not exphcitly appear in (P2). They are created immediately by the reflections at 
the walls and the piston at time r = 0. since the initial density 7r(?/, u,0) does not have 
to satisfy (H2)-(H3). Those discontinuity lines are y — 0, y — 0.5 and y = 1 at r = 0, 
and their images at r > will be slanted lines 

y — VT, y — 0.5 + vt, y — 1 + vt (2.34) 

respectively, see Fig. 2. So, their slope at any time r is positive and constant: dy/dv = r. 
It is not bounded away from zero as r — > 0, so we have a technical violation of (P2') for 
small r, but it will be clear immediately why this does not bother us. 

The singularity lines (2.34) only exist in the region Vmin < 1"^! < Vmax (elsewhere 
p = 0), hence they cannot intersect the piston y — Y{t) for small r. It will take some 
time, at least 

T* = — > 

'"max + l3eo 
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Figure 2: Slanted discontinuity lines. 
The dashed vertical line shows the piston position. 

before any of these singularity lines "reaches" the piston and its effect has to be reckoned 
with. At that time the slopes of those lines will be bounded away from zero: dy/dv > 
T* > 0, hence (P2') will hold. 

We now consider the evolution of all discontinuity curves of the function 7r{y,v,T) 
as T increases. Let a discontinuity curve of the function 7i{y,v,s) at time s be given by 
equation y = gs{v), and its slope is then hs{v) = dgs{v)/dv. Since the curve and its slope 
change in time, the function Qs and its derivative hg depend on s. According to (2.2), we 
have gs+r{v) = Qsiy) + vr between collisions with the piston and the walls, hence 

= , and ^^ = 1 (2.35) 
ds ds 

Hence, between collisions with the piston, the slope of discontinuity curves grows linearly 
with r (note that, in particular, it remains positive). 

Now, let the curve y = gs{v) cross the piston at some point 

9s{v) = Y{s) (2.36) 

(this equation makes v a function of s). After the collision with the piston, this point 
transforms to (Y,2W — v), according to the rule (H3), here W = W{s) is the piston 
velocity. If r > s is some fixed time, then the image of our point at time t is {Y + {t — 
s){2W — v),2W — v). Such points make a curve on the y,v plane, parameterized by s 
(the collision time). This will be the discontinuity curve for the density 7r(?/, v, r) at time 
r. Let ys = Y + {t — s){2W — v) and Vg = 2W — f be the coordinates of a point on that 
curve. To compute the slope dys/dvg of that curve, we first differentiate ys and Vg with 
respect to the parameter s: 



^ = W + {r-s) 
ds 



\^dW _ dv 
ds ds 



- {2W - v) 
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v-W 



2A 

ds 



and 



dvs ^dW dv ^ , dv 
= 2— -^2A- — 



ds ds ds ds 
where A = A{s) is the piston acceleration (at the colhsion time s). Also, differentiating 
(2.36) with respect to s and using (2.35) gives 



dgs{v) dv 



+ v^W 



hence 



dv ds 
dv W — V 



ds hs(v) 

Therefore, the slope of our singularity curve at time r is 

dy ^ {v - W) [hsjv) + T-s]+ 2Ahs{v){T - s) 
dv^^' V -W + 2Ahs{v) 

As long as (P1')-(P5') hold, we have W = O(£o) and A = O(£o), hence 

^(r) = /i,(^;)+T-s + 0(£o) 
dv 



(2.37) 



(2.38) 



Hence, every collision with the piston only adds a O^Eq) correction to the linear growth 
of the slopes of discontinuity curves. 

Next we check the conditions (P3')-(P5') based on the following lemma: 

Lemma 2.11 Let (P1')-(P5') hold on an interval (0,r). Then for every point {y,v) e 
Q~^{t) there is another point {yo,VQ) e Q'^ such that Tr{y,v,T) = 7i{yo,VQ,0) and 



l^o|| = 2«3,T+l)-B£o 



Proof. We set {yo,vo) = J-'^'^{y,v) and use (2.32). At each collision of the point {yo,vo) 
with the piston, its speed 1^1 changes by 2\W\ < 2B£q according to (2.31) and (2.26), 
and the number of collisions is bounded by v'^^-^t + 1. □ 

Lemma 2.11 immediately implies that (P3') and (P4') continue holding with global 
constants I'J, V2, v'^^^ and v'^^ slowly changing with time - they change at most by CTsq 
on the interval (0,T), with a global constant C > 0. In particular, v[ and v'^^^ remain 
positive, provided Sq is small enough. The constants ti'^i^ and tt^^x change at all. 

To check (P5'), we explicitly construct the map Rr '■ {y,v) {y*,v*) involved in 
(2.14) and (2.15), it is defined here by Rr = J-"^ o Rqo JF""", where Ro{y, v) = {1 — y, —v) 
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is a simple reflection "across the piston" at time r — 0. Now, (2.13) follows from (2.26), 
and (2.15) follows from (P5). 

Lastly, we derive (2.14) from the Lemma 2.11. Let {y,v) be a moving point at time 
r and {yo, vq) = J-'~'^{y, v) e its preimage to time zero. Compare the evolution of the 
point {yo,vo) and its mirror image i?o(|/o,^'o) = (1 ~ 2/o, — ^^o) G "across the piston" 
during the interval (0, r). Due to (2.13) and (2.26), these two points will experience 
collisions with the walls and the piston at time moments that differ at most by 0{eo). 
And their velocities will also differ at most by 0{eQ), hence their positions at time T will 
be almost symmetric about the piston, up to 0{so). This implies (2.14). 

Note that by the given time T the above two moving points may have experienced 
a different number of collisions, as one point may have just collided with the piston or 
a wall, while the other may be about to collide with it. To take care of this case, we 
introduced the special distance dy in (P5'). □ 

We summarize our main results in the following theorem: 

Theorem 2.12 Let T > be given. If the initial density TCo{y,v) satisfies (P1)-(P5) 
with a sufficiently small Eq, then 

(a) the solution of our hydrodynamical equations (H1)-(H4) exists and is unique on the 
interval (0, T); 

(b) the density 7T{y,v,T) satisfies (P1')-(P5') for all < t < T; 

(c) The piston velocity and acceleration remain small, \W{t)\ = 0{eo) and \A{t)\ = 
0{eo), and its position remains close to the midpoint 0.5 in the sense \Y{t) — 0.5| = 
0{so), for allO<T < T; 

(d) we have \Tk — k/v^axl — O{eo) for all 1 < k < Tv^ax, o-nd ifTvuiax > 2, then also 
\n - r2\ = 0{eo). 

Corollary 2.13 If Eo = 0, so that the initial density 7ro(y,f) is completely symmetric 
about the piston, the solution is trivial: Y{t) = 0.5 and W{t) = for all r > 0. 

Lastly, we demonstrate the reason for our assumption that all the discontinuity curves 
of the initial density 7r{y, v) must have positive slopes. It would be quite tempting to let 
n{y, v) have more general discontinuity lines, e.g. allow it be smooth for Vmin < 1^1 < Vmax 
and abruptly drop to at = Vmin and v — fmax- The following example shows why this 
is not acceptable. 

Example. Suppose the initial density Tro{y,v) has a horizontal discontinuity line v — Vq 
(say, Vq = finin Or Vq — I'max)- After ouc interaction with the piston the image of this 
discontinuity line can oscillate up and down, due to the fluctuations of the piston accel- 
eration (Fig. 3). As time goes on, this oscillating curve will "travel" to the wall and come 
back to the piston, experiencing some distortions on its way, caused by the differences 
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Figure 3: A horizontal discontinuity line (bottom) comes off the piston as an oscillating 
curve (top). 

in velocities of its points (Fig. 3). When this curve comes back to the piston again, it 
may well have "turning points" where its tangent line is vertical, or even contain vertical 
segments of positive length. This produces unwanted singularities or even discontinuities 
of the piston velocity and acceleration. The same phenomena can also occur when a 
discontinuity line of the initial density 7iQ{y,v) has a negative slope. 

3 Dynamics before the first recollision 

In this section we begin to study the mechanical model of the piston in the ideal gas 
described in Section 1. We will show that the random trajectory of the piston described 
by the functions Yl{t,u) = Xl{tL,u)/ L and Wl{t,u!) = Vl{tL,u!), cf. (1.6), converges 
in probability, as L ^ oo, to the solution of the hydrodynamical equations Y{t) and 
W{t) found in the previous section, on the zero-recollision interval (0,ri). 

Convention. For brevity of notation, we will suppress the dependence of L and u in 
our expressions, when it does not cause confusion. For example, we will write X{t) and 
V{t) instead of XL{t,uj) and VLit,u), respectively, etc. 

We will work here with the microscopic time t. First, we define the "microscopic" 
gas density, which we will denote by p{x, v, t), for all t > 0. For t = it is initialized by 
p{x,v,0) = Tro{x/L,v), see (1.7). For t > 0, its evolution is defined by the rules similar 
to (H1)-(H3): the free motion between collisions 

p{x,v,t) = p{x — vs,v,t — s) (3.1) 

for s > such that x — vr ^ {0, X{t — r), L} for all r G (0, s); reflections at the walls 

p{0,v,t) = p{0, —v,t) and p{L,v,t) = p{L, —v,t) (3.2) 
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and elastic collisions with the piston 

p{X{t)±Q,v,t) = p{X{t)±0,2V{t) -v,t) (3.3) 

Since the last equation involves the random functions X{t) and V{t), the density p{x, v, t) 
will depend on u, i.e. it is now a random function. 

The evolution of the density p{x, v, t) can be conveniently described with the help of 
a one-parameter family of transformations F* similar to J-"^ defined in Section 2. Let 
{x, v) be a point in the domain 

G := {{x,v) : < x < L} 

Its trajectory {xt, Vt) for i > is defined by the free motion inside the container, xt — Vt 
and Vt — whenever Xt ^ {0, X(t), L}, reflections at the walls {xf+o, Vt+o) — {xt-o, —Vt-o) 
whenever Xt_o G {0, L}, and collisions with the piston 

{Xt+O: Vt+o) = (Xt-O: '2V{t) - Vt-o) (3.4) 

whenever Xt-o = X{t). Now the family of transformations is defined by F^{xq,Vq) — 
{xt,Vt) for t > 0. We will also write F~^{xt,Vt) = {xq,vo). Now we simply have 

p{xt, Vt, t) = p{F'\xt, Vt), 0) = p{xo, vo, 0) (3.5) 

Note that for each t > the map F* is a bijection of G and preserves area, i.e. 

det\DF\x,v)\^l (3.6) 

We emphasize that the transformation F*, just as the density p{x,v,t), is random, i.e. 
depends on cu. 

Remcirk. The piston velocity V{t) is a piecewise constant function updated at the 
moments of coUision with gas atoms by the rules (1.2)-(1.3). If t is such a collision 
moment, then V{t) in equation (3.3) must be replaced by the average of its one-sided 
limit values {V{t — 0) -|- V{t + 0))/2. This modification is important, since it makes the 
rule (3.3) equivalent to (1.2)-(1.3) when {xt,Vt) represents an actual gas particle of mass 
m. Otherwise, it will correspond to the motion of a particle of zero mass, and we may 
call it a virtual particle. 

Because of (P4), the initial density p{x, v, 0) can only be positive in the region 

:= {(X, V) : <X < L, V^in < \v\ < Vmax} 

which therefore contains all the gas particles at time t — 0. For any t > 0, the region 
G~^{t) := F^{G^) contains all the actual gas particles at time t, and p{x, v,t) — Q for all 
{x,v)iG+{t). 



27 



For each point {x, v) & G and t > Owe define the number of coUisions with the piston 
during the interval {0,t) 

N{x,v,t) = #{s e (0,t) : Xs = X{s), + V{s)\ 

Then we partition the region G into subregions 

Gn{t):={F\x,v): {x,v)eG & N{x,v,t)=n} 

and put G^(t) := G~^(t) fl Gn(t). The region G^(t) is occupied by the points that have 
experienced exactly n collisions with the piston during the time interval (0,t). 

We emphasize that our transformations F* and the regions Gn{t) and depend 
on cu, i.e. are random. We note, however, that they are completely determined by the 
trajectory of the piston, i.e. by the function X{s), < s < t. 

The family of transformations J-"^ '■ Q ^ Q introduced in Section 2 induces a (de- 
terministic) family F* : G ^ G defined as follows: if J-"^{y,v) = {yTjVr), then we put 
F"^^{yL,v) := [y-rL^v). The transformations F* define a deterministic density function 
on G by p{x,v,t) :— p{F*{x,v),0), which is related to the density 7r{y,v,T) studied in 
Section 2 by ^(a;,^,^) = n{x/ L,v,t/ L). We also put G+{t) := J^{G+). 

The gas particles in Go(^) make a Poisson process, as the following lemma shows. Let 
uj & Vt and t > 0. Fix the trajectory of the piston X[s), < s < t. That completely 
specifies the region Go{t) and the density p{x,v,t). 

Lemma 3.1 The conditional distribution of the gas particles in Go{t) (given the tra- 
jectory X{s), < s < t, of the piston) is Poisson with density function L'^p{x,v,t). 

Proof. Let D C Go{t) be any domain. Then its preimage F~^*~^\D) stays positive 
distance away from the piston X{s) for all s G (0,t). Hence, the particles starting out in 
the region F~^{D) and ending up in the region D could not affect the piston during the 
time interval (0, t). Therefore, the number of particles in D at time t, being equal to the 
number of particles in F~^{D) at time 0, is independent of the piston trajectory, so it is 
a Poisson random variable with parameter 

The identity of the above integrals follows from (3.5) and (3.6). □ 

Remark. For any domain D C Go{t) its preimage F~^{D) is actually independent of oo. 
Indeed, let Fq be another family of transformations on G defined by the free motion on 
the entire interval < x < L and elastic reflections at the walls x — and x — L only (as 
if the piston did not exist). Then we have F''*{D) — Fq*{D) for any domain D C Go{t). 
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For n > 1, we define Tn to be the earfiest time the piston interacts with points from 
(thus creating the region G^+i{t)), or, equivalently, 

7; = sup{G++,(t)=0} (3.7) 

The time moments T„ = Tk(a;) are random analogues of t„ introduced in Section 2. In 
particular, Ti is the time of the first recoUision in the system (by an actual or a virtual 
particle) . 

Lemma 3.2 For all uj & Q 



Proof. The fastest particles {x, v) G that collide with the piston at time will move 
with the speed fmax and recoUide with the piston at time t that satisfies Wmax^ + l^(^) ~ 
L/2\ — L. For all such t we have 

^max 

This proves the lemma. □ 

During the time interval (0, Ti) the piston interacts with particles in Gq (t). Denote 

by 

^o{t) = {{X, v) : X{t) + 0, -fmax <V < -Vmin} 

U{(X, V) : X^ X{t) - 0, Vrnin <V < V^a^} 

two immediate one-sided vicinities of the piston which contain all "incoming" particles, 
which are about to collide with the piston. 

In order to study the piston dynamics on the zero-recoUision interval, it is convenient 

to assume that the piston is slow enough and only interacts with the original particles 
that started in at time 0. A subinterval (0, 5*1) C (0, Ti) where the piston satisfies 
these requirements will be called a "slow" interval: 

Definition We define (0, 5*1) C (0, Ti) to be the maximal time interval on which 

(a) \V{t)\ < 

(b) Mt) c Go{t). 

Note that by the condition (b) all the incoming particles that are about to hit the 
piston at time t have started out in at time zero and have never interacted with the 
piston during the time interval (0, t). 

Almost all of our considerations in this section are restricted to the "slow" interval 
(0, Si). But in the end of the section we will see that for typical uu & Q the "slow" interval 
coincides with the entire interval (0,Ti), i.e. = Ti. 
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Now, for every t e (0, Si) we define the density of colliding particles "on the piston" 
q(v, t) by 

= I "'J! (3.9) 

' ^ [ p{X{t) + 0,v,t) if v < ^ ^ 

cf. (2.17). Next, we define 

Qo{t) = J {sgn v)q{v,t)dv (3.10) 

Qi{t) = J v{sgiiv)q{v,t)dv (3.11) 

Qiit) ^ J v'^{sgnv)q{v,t)dv (3.12) 

in a way similar to (2.19)-(2.21) in Section 2. 

Since p(x, v, t), restricted to the domain Go(^), coincides withp(a;, v, t) — 7r{x/L, v, t/L), 
the conditions (P1')-(P2') imply 

Lemma 3.3 The density p{x,v,t) restricted to the region Go{t) is piecewise smooth 
on the x,v plane. We also have \dp{x,v,t)/dx\ < D[/L. The discontinuity lines of 
p{x, V, t) within the region Go{t) have slope of order 0{1/L) (for large L, they are almost 
parallel to the x axis). 

Due to the above lemma the quantities Qj, i = 0, 1, 2, are, as functions of the piston 
position X, smooth and have derivatives 



^Q^ 



dX 



< ^ (3.13) 



where const is a global constant. 

The following theorem gives a key technical estimate of this section. 

Theorem 3.4 For sufficiently large L there is a set 0,^ (Z ^ of initial configurations of 
gas particles such that 

(i) for some constant c > 

P(fi*) > 1 -L-'^i'ii"^' 

(ii) for each configuration a; e Oq, for each time interval 

{t,t + At) C {0,Si) 

such that 

the change of the velocity of the piston during {t, t + At) satisfies 

V{t + At) -V{t) = V{t)At + x (3.15) 
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where 

V{t)/a = Qo{t)V^{t) - 2Qi{t)V{t) + QsW (3.16) 

and 

InL VAt 

IXl < ^' 

with some constant C > 0. 



X\<C ^ (3.17) 



Remark. The function V(t)/a in (3.16) is the random analogue of the quadratic polyno- 
mial (2.18). The term V{t) At in (3.15) is the main ("deterministic") component of the 
dynamics of the piston velocity. The term x represents random fluctuations. 

Proof. The set VIq will consist of all configurations that satisfy certain requirements. We 
start with preliminary requirements. 

Consider a discrete set of time moments ti — i/L?', where i — 0, 1, . . . , / and / = 
[7i,max-t'^]. Partition the domain G into the strips Sj := {{x,v) : j/L'^ < x < (j + l)/L^}, 
where j = 0, 1, . . . , — 1. For each i and j denote by Nij the number of gas particles 
in the region Sj fl Go(^i) at time ti. Our preliminary requirements are 

Nij<\nL (3.18) 

for all < i < / and < j < L^. We observe that Nij equals the number of gas particles 
in the region F~**(Sj fl Go{ti)) at time 0, and 

F-'^iS,nGoiti))cF,''{S,) 

So, Nij does not exceed the number of gas particles in Fq^^^Sj) at time 0, denote the 
latter by N^j. Now, Nij is a Poisson random variable whose parameter Xij is bounded 

by ' ' ^ 

— ^maxI'S'jl ^ 27rmax(^max ^min) 

According to Corollary A. 4, for each i,j our requirement (3.18) will fail with probability 

< /.-f^ininL -^j^]^ some d > 0. The total number of pairs i,j equals L^I = L^Ti ^ 
L^/2fmax. Hence, all our preliminary requirements hold with probability > 1 — L~'^ ininL 
with some global constant c' > 0. 

We now turn to the proof of (3.14)-(3.16). Let i = [LH] and ti = (i + l)/L^. Note 
that ti — t < L~^. One can easily derive from our preliminary requirements that the 
number of gas particles colliding with the piston on the time interval {t,ti) is less than 
const- In L. Hence, the piston velocity V does not change by more than const- In L/L^ 
during this interval. This amount is less than the bound on x in (3.17) for all At satisfying 
(3.14). Therefore, we can ignore the interval (t, ti) and assume that t = ti. Note that the 
quantities Qo, Qii Q2 in (3.16) will change, as the result of the substitution t — ti, also, 
but only by the amount < const/L^ due to Lemma 3.3. This change can be obviously 
ignored, too. So, we suppose that t = i/L'^ for some i = 0,1, . . . , I 

Now we state our main requirements. We again partition the x axis into intervals 
Sj := {j / L"^ <x<{j + 1)/L^}, where j = 0, 1, . . . , — 1. For each j we put Xj = j/L"^. 
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For each integer p, \p\ < VuiinL^, we put Vp — p/ LF' and for each integer 1 < g < we 
put dq = q/ LP'. For each triple q) we now define two trapezoid-hke domains on the 
x,v plane (see Fig. 4 below): 



^'^'^ I X- Xj dq 



<v < v^^ } (3.19) 



and 



D^pn ■■= I : ^^-^ < -v^^ <v< -v^A (3.20) 

■''^'^ [ X — Xj dq J 

(here —l/dq is the slope of the oblique side of these trapezoids), and two strips 

U~p,q {{X, V) : \X - Xj + dq{v - Vp)\ < lOv^a^/L^, Umin <V< fmax} 

and 

Uj^p^q ■= {{X, V) : \X- Xj + dq{v -Vp)\ < lOv^a^/L^, -Vmax <V < -fmin} 

Note that Uf^p^q are the neighborhoods of the oblique sides of the trapezoids Dfp q. 

Consider all time moments ti = i/L'^ for i = 0,1, . . . , I. Denote by N^j^p^q the number 
of particles in the region FQ^'{Dfp q) at time 0. And denote by M^j^p^g the number of 
particles in the region FQ^'-iUfpq) at time 0. These are Poisson random variables. The 
parameter of the variable N^j^p^q is 

>^tp,,-E{NM,p,,)-L'J pix,v,0)dxdv 

One can verify directly that )<t,j,p,q uniformly bounded below by a positive constant 
(even for the smallest dp, i.e. for dq = l/L"^), due to the assumption (1.8) on the initial 
density. Also note that the parameters of are uniformly bounded above, by 

Our main requirements are 

I^4p..-^£,J<1^^\AS^ (3-21) 

and 

M4,,^<lnL (3.22) 

By Lemma A. 3 and Corollary A. 4 the probability that any of these requirements fails 
will be less than i^-<^inin^ with some constant d > 0. The total number of quadruples 
{hj:P:Q) does not exceed -L^Ti^max < L"^" with some fixed A" > 0. Therefore, all our 
main requirements hold with probability > 1 — L~'^ ininL ^j^jj some constant c" > 0. 
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In addition, let 

taken at time 0. This is an "integrated" Poisson random variable, as defined in Appendix. 
(Technically, we require there that the domain must be on one side of the x axis, and 
now it may happen here that the region Fq~** {Dfp g) crosses the wall x = or x = L and 
then lies on both sides on the x-axis; in that case we need to replace ti by a nearby time 
moment tj/ < ti so that Fq*^' {Df^,^ lies entirely on one side of the x axis and define 
^t.j;p,q time ti — til rather than 0; Some obvious modifications need to be made then, 
we omit details.) The estimates obtained in Appendix yield 



EiZtj,J = L'J vp{x,v,0)dxdv (3.23) 

j,p,q' 

Var(Z±.p^^) = L2 / ^ vMx,v,0)dxdv (3.24) 
Our last main requirement is 



- ^(Zj,.,J| < InL V^Var(Z±.,J (3.25) 

for all i,j,p,q. The probability of failure for these requirements is estimated exactly as 
above, by using Lemma A. 6. 

We now turn to the estimation of the piston velocity on the time interval {t,t + At). 
Recall that t — ti for some i = 0, 1, . . . , /. 

Velocity decomposition scheme. Here we obtain a general formula for the piston 
velocity, which we will use in the proof of several theorems. The laws of elastic collisions 
imply [LPS] 

k 

V{t + At) = (1 - e)''V{t) + £ E(l - ^)'"' ■ (3-26) 

i=i 

Here k is the number of particles colliding with the piston during the time interval 
{t,t + At)j and vj are their velocities numbered in the order in which the particles collide. 
Equation (3.26) can be easily verified by induction on k. 
We modify the formula (3.26) as follows: 

k 

V{t + At) = (1 - ek)V{t) +ef2vj + X^'^ + X^'^ (3-27) 

where 

X^^) = Vit)[{l-e)'-l+ek] 

and 

X^'^-ej:vA{l-e)'-^-l] 
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Let us assume that the fluctuations of the velocity V{s) on the interval {t,t + At), are 
bounded: 

sup \V{s)-V{t)\<6V (3.28) 

Consider two regions on the x, v plane: 

f/ X V -V(t) - (sgnv)6V 1 , , 1 , , 

Dl = {X, V) : ^ ' ^ < - — , Vrain < \v\ < t^max (3.29) 



x~X{t) At'' 



and 



f/ X V -V(t) + (sgnv)SV 1 , , 1 , , 

D2 = Ux,v) : j-xit) ^ "At' ^ ^ """^""J ^ ' ^ 

Each of them is a union of two trapezoids Di = Df U , i = 1,2, where denotes 
the upper and Df the lower trapezoid, see Fig. 4. 




m 

X 




Figure 4: Region Di is bounded by solid lines. Region D2 is bounded by a dashed line. 

Note that Di C -D2- The bound (3.28) implies that all the particles in the region 
Dl necessarily collide with the piston during the time interval (t,t + At). Moreover, the 
trajectory of every point {x,v) G Di hits the piston within time At, hence Di C Go{t) 
by the condition (b) in the definition of the "slow" interval (0,5*1). The bound (3.28) 
also implies that all the particles actually colliding with the piston during the interval 
(t, t + At) are contained in D2 (though it is not necessarily true that D2 C Go{t)). 
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We now obtain an upper bound on k. Since the velocities of the gas particles and the 
piston arc bounded by ^Jmax all the particles colliding with the piston during the interval 
(t, t + Ai) are contained in the region S' fl Go(^), where 

S'^{{x,y): \x-X{t)\<2v^^M} 

Hence, k does not exceed the number of particles in the region Fq~*(S") at time zero. Our 
preliminary requirements imply 

k < 2vinaxAt In L (3.31) 

This allows us to bound the quantities x^^^ ^iid X^^^ i for large L, by 

IX^^^I < v^^^e'k^ < Aa'vl^ilnLf (At)^ (3.32) 

for w = 1, 2. It also follows from (3.27) that 

SV < baVme^^At In L (3.33) 

We denote by k^ the number of particles in the regions n Go{t) for r = 1, 2 at time 
t = ti. We also denote by k~ the number of particles colliding with the piston "on the 
left", and by k'^ that number "on the right" (of course, k~ + k'^ — k). Clearly, 

ki <k^ < kf (3.34) 

Since Di C GQ{t), then kf equals the number of particles in the region FQ^(Df) at 
time 0. Similarly, kf equals the number of particles in the region F^^^Df n Go{t)) at 
time 0. 

The trapezoids Df, r — 1,2, can be well approximated by some trapezoids Dfp^^ 
defined earlier in our main requirements. In fact, the horizontal sides v — ii^min and v — 
if max are common for all trapezoids, only the vertical side x = X(t) and the oblique side 
need approximation. The symmetric difference between Df and the approximating Df^p^^ 
will lie inside some strips Sj and U^p^g also defined above. Then our main requirements 
will guarantee that 

kf > Xf-lnlJ)^ 



and 
where 

and 



kf < + lnL^A^ 
Xf = L^f p(x,v,0)d.dv=L^f p(.,v,t)±cdv 

JF-\Df) JDf 

\f — pix, V, 0) dx dv 

JF-*(D±nGoW) 

< Xf + L^7r^,^\Df\Df\ 
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where we used the boundedness of the density (1.9). The area of the domain D2 \ Df is 
bounded by 

\Dt\Df\<Av^^8VM 

Also, (1.9) imphes that 

Combining the above estimates gives 

|A;^-A^| < const -{LlnLVAt + L^SV At) (3.35) 
We now turn to the quantity 

k 

^ = E 

also involved in the main equation (3.27). Again we can decompose Z = Z~ + Z^, where 
Z~ and Z^ denote the sum of velocities of the particles colliding with the piston "on the 
left" and "on the right" , respectively. We put 

{x,v)eD^ 

(taken at time t). Analysis similar to the previous one and the requirement (3.25) with 
formulas (3.23)-(3.24) give 

\Z^-E{Z^)\ < const ■{L\nLVAi+L'^5V At) (3.36) 

We now combine (3.27) with all the subsequent estimates and obtain 

V{t + At) - V(t) = -sXiV{t) + eE{Zi) + x^^^ (3.37) 

where 

Xi — Xf + X]^ = / p(x, V, t) dx dv 

JDi 

E(Zi) = E(Z^) + E{Z^) = L'^ I vp(x,v,t)dxdv 

JDi 

and 

Ix^^^ I < const ■ [l-^ In L \fAi + 5V At + (In L)^ (At)^" 

Using Lemma 3.3 (as we already did in deriving (3.13)) to bound possible fluctuations 
of the density p(a;, t) within Di gives 

/ p{x,v,t)dxdv^{Qi{t)-QQ{t))V{t)At + x^^^ (3.38) 
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and 

/ vp{x,v,t)dxdv^{Q2{t)-Qi{t))V{t)At + x^^^ (3.39) 

with 

IX^"^ I < const • [(At) VL + 5V At] (3.40) 

for u — 4:,5. 

Therefore, we get 

V{t + At) - V{t) = V{t) At + x (3.41) 

where 

Ixl < const • [l-^ In L y/At + SVAt+ (In Lf {Atf] (3.42) 

By using (3.14) and (3.33) it is easy to incorporate the second and the third terms in 
(3.42) into the first one. Theorem 3.4 is proved. □ 

For the next theorem, we rewrite (P5) in terms of the microscopic coordinates: 

\p{x,v,0)-p{L-x,-v,0)\<eo (3.43) 

We also assume also that the initial velocity of the piston is set to V{0) — W{0) rather 
than to 0, see a remark Section 2 and another remark below. 

Theorem 3.5 Assume that Eq > is small enough. For all sufficiently large L, for each 
configuration cu & ft^ and for all t e (0, Si) we have 

(i) there is a constant B > such that 

\V{t)\<Beo (3.44) 

(ii) there is a constant Cq > such that 

mt) - V,(t)\ < ^°<^ (3.46) 

where Vo{t) is defined by 



Vo{t) = ^ . (3.46) 



whenever Qo{t) ^ and by 



otherwise. 
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Proof. We will start proving (i) and obtain (ii) as a "side result" . 

Assume that (i) is false and let < Si be the first time (3.44) fails, i.e. let 

\V{U + 0)\>B6o (3.48) 

Since (3.44) holds for t < t^, the piston's position satisfies 

\X{t) - L/2\ < Beot (3.49) 

for all t < t^. Assume for the moment that the piston did not move at all, i.e. X{s) = L/2 
for all < s < t. Then, by (3.43), the density q{v, t) of the gas "on the piston" would be 
almost symmetric, i.e. 

\q(v,t) - q(-v,t)\ < £o 

Hence, we would have 

\Qi{t)\ < CiEo, i = 0,2 (3.50) 

with Ci = ■I'max- We emphasize that Ci does not depend on the choice of B in (3.44). 
Below we will introduce some more constants C,, i > 2, so that none of them will depend 
on B. 

When the piston actually moves and covers the distance X{t) — L/2, then (3.13) and 
(3.50) imply 

|Q,(t)| < Ci£o + C2L-i|X(t) -L/2|, i = 0,2 (3.51) 

with some constant C2 > 0. At the same time, Qi{t) stays bounded above and below by 
positive constants for all t < t^: 

< (9l,min < Ql{t) < (3l,max < OO (3.52) 

where (5i,min and (5i,max are constants determined by TTmin and TTmax in (P3)- 
By (3.49), (3.51) and (3.8) we have 

\Qi{t)\<{Ci + Bv-l,)eo, t = 0,2 (3.53) 

for all t < t^. This and (3.52) imply that the quadratic polynomial (3.16) has two real 
roots, and the smaller one (the one closer to zero) is given by (3.46)-(3.47), cf. (2.24) 
and (2.25) in Section 2. Due to (2.22) and (3.52) we have 

\Vo{t)\<C-iSo + C^L-'\X{t)-L/2\ (3.54) 

for all t < and some constants C3, C4 > (independent of the choice of B). 
We note that (3.54) and (3.8) imply the boundedness of Vo{t) 

\Vo{t) I < const = CsEo + C^v-l, (3.55) 

for all t < t^. 
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Next, we need to estimate the derivative dVQ{t)/dt. (Note: the function Vo(t) defined 
by (3.46)-(3.47) is continuous and piecewise diffcrcntiable, it should not be confused with 
the piston velocity V{t), which is piecewise constant and hence not even continuous). Due 
to (3.13) 



dQi{t) 



dt 



dQi{t) 



dX 



V{t) 



< 



const • £o 



(3.56) 



for i = 0, 1, 2. Differentiating the quadratic equation 

Qo{t)Vo\t) - 2Q,{t)Vo{t) + Q2{t) = 

with respect to t gives 

dVo{t) ^ {dQo/dt)Vi - 2{dQi/dt)Vo + {dQ^/dt) 
dt ~ 2(gi - QoVo) 

Due to (2.22), (3.52), (3.55) and (3.56) we have 



(3.57) 



dVo{t) 



dt 



< 



(3.58) 



where > is a constant. 

Now consider the quantity T> defined by (3.16) as a function of V (with fixed Qi, 
i = 0, 1, 2). Its derivative is 

dV 

— = 2a[QoV-Q,] 

Due to (3.52) and (3.53) there are positive constants < Ei < E2 such that for all t <t^ 
we have 

^ dV 

-E2<^< -E, 

therefore, by the mean value theorem, for all t <t^ 

V{t) 



-Eo < 



Vit) - Vo{t) 



< -El 



(3.59) 



We now prove (3.45) for all t < t* with some constant Co > independent of the choice 
of B in (3.44). Recall that wc have set the initial velocity of the piston to V{0) — W{0) 
and that 14^(0) = Vo{0), see Section 2. Now, by way of contradiction, let t e (0, i*) be 
the first time (3.45) fails. Denote by 



1 



" L2/3 1nL 

the maximal allowed time increment in Theorem 3.4. Let s — t — Ag. Due to Theorem 3.4 

V{t) = V{s)+V{s)Ao + x (3.60) 
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with 

C^/A^\nL _ C^ML 

IXl < I - 2.4/3 

Due to (3.58) we have 

Voit) = Vois)+xo (3.61) 

with 

IXol < 



For brevity, put U{s) — V{s) — Vo{s) for all s. Subtracting (3.61) from (3.60) then gives 

U{t) ^U{s)+V{s) Ao + x' (3.62) 
with x' = X ~ XO) so that for large L 

1^ I ^ ^4/3 (3-63) 

Now assume, without loss of generality, that U{t) > 0. Since (3.45) fails at time t, we 
have 

Ui*) >. ^ (3.64) 
Now consider two cases. If U{s) < 0, then by (3.59) 

fin L)^/^ 

Uit) < \Vis)\Ao + \x'\ < E2 ms)\Ao + \x'\ « ^ 
for large L, which contradicts to (3.64). If U{s) > 0, then, again due to (3.62) and (3.59), 

U{t)<U{s)[l-E^Ao]+x', 

hence 

^^'^ ^ rr^>(^w-^')(i+^iAo) 

> U{t) + U{t) El Ao-2x' (3.65) 

Now, if Co in (3.45) is large enough, say Cq = 5C/Ei, then U{f)EiAo > 2x' by (3.64) 
and (3.63). This fact and (3.65) imply U{s) > U{t), so (3.45) fails at an earlier time 
s < t, a, contradiction. Hence, (3.45) is proved for all t < and Co = 5C/Ei. 
Now, combining (3.54) and (3.45) gives, for large L and all t < U 

\dX{t)/dt\ < ICsEo + C4L-^\X{t) - L/2\ 

Using the standard Gronwall inequality in differential equations, see, e.g.. Lemma 2.1 in 
[TVS], gives 

\X{t) -L/2\< 2£oC3C4-iL(e^*^"'* - 1) 
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and 

\V{t)\<2eoCse^''^-'' 

for all t < t^. By (3.8) we have 



< 2£oC3e^^^-^ (3.66) 

Now we choose B — 3C3e^*^™^. Clearly, (3.66) then contradicts (3.48). This com- 
pletes the proof of (3.44) for all t < Si. Theorem 3.5 is now proved. □ 



Corollary 3.6 Assume that Eq > in (1.11) is small enough. Then, for all large L and 
all uj G l^Q, we have S\ = Ti, i.e. the previous theorems hold on the entire zero-recollision 
interval (0, Ti). 

Proof. Recall that the "slow" interval (0, 5*1) C (0,Ti) is defined by two conditions, (a) 
and (b). If 5*1 < Ti, then either (a) or (b) fails at Si. Clearly, (a) cannot fail "abruptly" 
since (3.44) holds for all t < Si. 

Suppose (b) fails at some s < Ti, while (a) still holds. The failure of (b) means 
that at time s the piston "collides" with a point {x,v) such that fmin < I'J^I < '^^max and 
{x,v) ^ Go{s). Therefore, the backward trajectory F~^'^~''^\x,v), < t < s, of the point 
{x,y) hits the piston at some time t > 0. Now, during the time interval (0,t) the piston 
covers the distance \X{t) — L/2\ < Bsot, and during the time interval (t, s) the trajectory 
of our point covers the distance |v|(s — t) < Vmax(5 — t). Hence we have 

L < \X{t)-L/2\ + \v\{s-t) + \X{s)-L/2\ 

< \X{s) - L/2\ + Vma^S - (t^max - Beo)t 

On the other hand, since s < 71, we have 

L> \X{s)-L/2\ + 

see the proof of Lemma 3.2. This contradiction shows that (b) cannot fail either. The 
proof of Corollary 3.6 is completed. □ 

Remcirk. We have reset the initial velocity of the piston to V{0) — W{0) here, while 
in Section 1 it was set to zero. If 1^(0) = 0, then Theorem 3.4 would imply that 
V(t) converges to Vo{t) exponentially fast in t, until it gets 5-close to Vo(t) with S = 
Co(lnL)^/^L~^/^. After that all our results will apply without change. The initial inter- 
val on which the convergence takes place will be of order In L, and in the hydrodynamical 
time it is L"^ InL, which vanishes as L — > 00. This is why we simply opted for the most 
convenient setting V{0) — W{0) here. 

The following theorem improves the results of Theorems 3.4 and 3.5. 
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Theorem 3.7 Assume that Eq > in (1-11) is small enough. Then there is a constant 
C > such that for all large L, all a; e and all t < Ti 

\V{t)-Vo{t)\<^^ (3.67) 
and for any time interval {t,t + At) C (0,Ti) such that < At < 1 we have 

\V{t + At)-V{t)\ ^c ^^^^f^ (3.68) 



Proof Due to our choice of the initial velocity, V{0) = Vo(0), hence (3.67) holds for at 
least small t. Assume that it fails at some < 71, and t* is the earhest time of failure. 
Without loss of generality, assume 

V{t,) - Vo{t,) > (3.69) 
Let < to < t^ he the latest time when 

^(^o) - Vo{to) < (3.70) 

Then we have 

Vito) - Vo{to) < V{t) - Vo{t) < V{t,) - Vo{t,) (3.71) 
for all t e (to, t*)- Let 

At = min{l,t -to} (3.72) 

Wc will analyze the dynamics of the piston during the time interval (to, to + At). Due to 
(3.58) we have 

\Vo{t) - Vo{to)\ < 5V := EoEoL-'At (3.73) 
for all t G (to, to + At), hence (3.71) implies 

V(t) > V{to) - 5V (3.74) 

We note that At is not too small, it is at least At > (L^/^lnL)~^. Indeed, otherwise we 
would have t* = to + At and then (3.69), (3.70) and (3.73) would imply 

V{t^) - V{to) > 2-^CL-HnL- EqEoL-^'^ 

which would contradict Theorem 3.4, since X>(to) < (because ^(to) > V^)(to)). 

Next, we develop a generalized version of the velocity decomposition (3.27) in the 
proof of Theorem 3.4. We partition the interval (to, to + At) into subintervals of length 
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5 (to be chosen shortly) with endpoints U — to + i5, i — 0,1, I , where / = At/S. We 
select S so that 

and At/ 5 is an integer (for convenience) . Preliminary requirements in the proof of The- 
orem 3.4 allow us to adjust time so that L^to L'^^ are integers, hence L'^ti will be an 
integer for every i. The velocity decomposition in the proof of Theorem 3.4 now applies 
to each subinterval (ij,ti+i) of length 5. In particular, (3.27) implies 

k 

V{ti+,) - V{ti) = -skiV{ti) +ef2vj + X?^ (3.76) 

where ki is the number of particles colliding with the piston during the time interval 
(tj,ti+i) and Vj, 1 < j < ki, are their velocities. The fluctuation term Xi^^ can be 
bounded by (3.32): 

\x^^\<8a\l^{\nLr5' (3.77) 
Due to (3.74), the expansion (3.76) can be rewritten as 

k- 

V{ti+,) - V{ti) < -ehV{to) +eJ2vj + xf ^ + xf (3-78) 

with 

Ixf ^ I < ski < 2aVraa^Eo EqL-^ In L At 5 
where in the last step we used (3.31) and (3.73). Summing (3.78) up over i yields 

k 

V{to + At) - V{to) < -ekV{to) + £ E + X^'^ (3-79) 

i=i 

where k is the number of particles colliding with the piston during the time interval 
{to, to + At) and Vj, I < j < k, are their velocities, and we have 

IX^'^I < 8a''vl,^{\nL)''6At + 2av^^EoeoL-HnL{At)^ 

< {Sa''vl^ + 2aVraa.Eoeo)L-HnLAt (3.80) 

(where we used (3.73) and the assumption At < 1). 

The expansion (3.79) can be analyzed similarly to (3.27) in the proof of Theorem 3.4. 
Define a region on the x, v plane: 

f, , v-V(to) + 5V 1 , , ] , , 

Di = <{X,V) : ^ _ j^^^^^ < --^^ ^min < \V\ < Vraa^ | (3.81) 

It is the union of two trapezoids Di = Df U Di , where Z^f denotes the upper and Df 
the lower one, see Fig. 4. The bound (3.74) implies that all the particles in the region 
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necessarily collide with the piston during the time interval {to, tg + At) and all the 
particles actually colliding with the piston on its left hand side during this interval of 
time are contained in Di . 

Since v > V{to) for all particles {x,v) G and v < V{to) for all {x,v) e Df, 
we can remove from (3.79) the particles that do not belong in and simultaneously 
add to (3.79) the particles that belong in but do not collide with the piston. This 
modification only makes the right hand side of (3.79) larger, hence 

V{to + At) - V{to) < -ehV{to) + sJ2'"j+ X^'^ (3-82) 

where ki is the number of particles in Di at time to, and the summation runs over all 
those particles. Let Zi = J2{x,v)€Di ^- Just like in the proof of Theorem 3.4, our main 
requirements stated there guarantee that 

- E{ki)\ < csLlnlVAt 

and 

|Zi - E{Zi)\ < CiLl-alVAt 

where the constant 03,04 > do not depend on the choice of C in (3.67), which we 
have not made yet. Now, computing the mean values of ki and Zi as in the proof of 
Theorem 3.4 we arrive at 

V{to + At) - V{to) < V{to) At + x^') + X^'^ (3.83) 

with 

IX^^^ I < «(C3 + C4)L-^ In L \fAi + C5 5V At (3.84) 

where C5 > is a constant independent of the choice of C in (3.67). The last term in 
(3.84) comes from the adjustment 5V to the velocity V{to) in (3.81). This last term is 
bounded by Cr^EQeoL'^^At)"^ , and since At < 1, it can be incorporated into the first term 
in (3.84). Recall that ^(to) — K)(^o) ~ 'i~^CL~^ In L > (here we have an approximation 
up to a quantity of order since the piston velocity changes by 0(1/L^) at each 

collision). Then due to (3.59) we have 

V{to) < -Ei{V{to) - Vo{to)) ^ -2-^CEiL-HnL 
Therefore, combining the above estimates gives 

V(to + At) - y (to) < -2-^CEiL-^ In L At + x^^^ (3.85) 
with x^^^ = X^^^ + X^^^ bounded by (3.80) and (3.84): 

Ix^^^l <C6{L-HnLAt + L-^\nL^/At) 
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where ce > is a constant independent of the choice of C in (3.67). Now we chose the 
constant C there as 

C = max{c6, QE^^cq) (3.86) 

Then (3.85) imphes 

V{to + At) - V(to) < cqL-^ In L (-2 At + 

and hence, due to (3.73), 

V(to + At)-Vo(to + At) < V(to)-Vo(to) 

+C6L-^ In L {-2At + VAt) + Eq e^L-^At (3.87) 

We now have two cases. First, let t^—to < 1, hence At < 1. The expression —2At + \/Ai 
has a maximum, equal to 1/8, at the point At — 1/16. Therefore, (3.87) imphes 

V{t,) - Voit,) < V{to) - Voito) + 8-%L-' In L + Eo £o^"' 

This contradicts (3.69) and (3.70) when L is large enough, recall our choice of C in (3.86). 
Consider the second case: t^ — to > 1. Then At = 1 and (3.87) implies, for large L, 

V{to + 1) - Vo{to + 1) < V{to) - Vo{to) - 2-^ceL-HnL 

which contradicts (3.71). This completes the proof of (3.67). 

We now prove (3.68). If At < (L^/^lnL)~^, then we can use Theorem 3.4: 

\V{t + At) - V{t)\ < \V{t)\At + CL-HnLy/Ai 
The early estimates (3.59) and (3.67) imply 

\V{t)\ < E2\V{t) - Vo{t)\ < CE2L-HnL (3.88) 

so that (3.68) follows (with some larger value of C than above). 

Now let (L^/^lnL)"^ < At < 1. Without loss of generality, assume that V{t + At) > 
V{t). Moreover, we can assume that 

V{s) > V{t) for all se{t,t + At) (3.89) 

Indeed, if this is not the case, we can replace t hy t' ~ max{s < t + At : V{s) < V{t)} 
and prove (3.68) for the smaller interval {t',t + At). 

Next, our plan is to apply some estimates from the proof of (3.67) and then argue along 
the lines of the proof of Theorem 3.4. Denote to = t and partition the interval (to, to + At) 
into subintervals of length 5 satisfying (3.75). Then we again have decomposition (3.76)- 
(3.77). Due to (3.89) we have V{ti) > V{to) for all i, hence (3.76) implies 

y (ti+i) - V{ti) < -ehV{to) +eY,Vj + xf^ 
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Summing this up over i gives (3.79) with the bound (3.80), in which the second term 

(2) 

can be simply removed, since we do not have xl anymore. Next, possible fluctuations 
of the piston velocity V{s) during the time interval {to, to + At) can be estimated with 
the help of (3.67) and (3.73): 

\V{s) - V(to)\ < 5V := 2CL-^ In L + Eq SoL-^At 

for all s e {tojto + At). Then we estimate the random variables k and Z = J2j "^j along 
the lines of the proof of Theorem 3.4, starting with construction of two domains Di and 
D2 by (3.29)-(3.30), etc. Repeating the argument almost word by word we arrive at an 
analogue of (3.41): 

V{t + At) - V{t) < V{t) At + x' 

where 

Ix'l < const • [l-^ In L VAt + SV At + L'^ In L At 

where the last term comes from (3.80), which we have now, instead of (3.32) (note: (3.32) 
would not be nearly enough anymore, since At is large; this is why we needed to partition 
the interval {t,t + At) into smaller subintcrvals). We now combine the above estimates 
with (3.88) and complete the proof of (3.68) and Theorem 3.7. □ 

We finally prove the convergence, as L — > cxd, of the random trajectory of the piston 
to the solution Y{t), W{t) of the hydrodynamical equations found in Section 2. 

Theorem 3.8 Assume that Eq > in (1.11) is small enough. Then, for all large L and 
all LU e Qq, there is a constant C > such that 

\YLiT,uj)-Y{T)\<^^ (3.90) 

and 

\WL(r,u;)-W(r)\<^-^ (3.91) 



for all < T < min{ri,Ti/L} and 



\n-T,/L\<^ (3.92) 



Proof In Section 2 we defined the function F{Y, r) so that the hydrodynamical solution 
Y{t) satisfies 

dY{T)/dT = F{Y,t), Y{0) = 1/2 (3.93) 
see (2.29). Now Theorem 3.7 implies that for all a; e JIq random trajectory satisfies 

dYL{T,u)/dT = F(y,T) + x(r,a;), ^^(0,0;) = 1/2 (3.94) 
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with some 



X{r,uj)\ < 



C InL 
L 



Recall that \dF{Y,T)/dY\ < k, see (2.30). Therefore, the difference Zl{t,uj) := Yl{t,lj)- 
Y{t) satisfies 



and Zl{0,uj) = 0. By the standard Gronwall inequality in differential equations, see, 
e.g.. Lemma 2.1 in [TVS], we have 



for all T < min{ri, Ti/L}, which imply (3.90) and (3.91). 

Lastly, we verify (3.92). By (3.91), random fluctuations of the piston velocity are 
bounded by CL~^ InL. Hence, random fluctuations of the velocities of particles that have 
had one collision with the piston are bounded by 2CL~^ In L. The random fluctuations of 
the positions of both the piston and particles at every moment of time t < min{riL, Ti} 
are bounded by the same quantities (with, possibly, a different value of C) in the coor- 
dinate y — x/L. On the other hand, the relative velocity of the piston and the particles 
stays bounded away from zero (by, say, fmin — ABeq > 0). Hence the time of the flrst 
recoUision Ti/L can differ from ti by at most const-L~^ InL. Theorem 3.8 is proved. □ 

4 Dynamics between the first and second recollisions 

In this section we study the one-recoUision interval (ri,r2), on which gas particles expe- 
rience the second collision (i.e., the first recoUision) with the piston. 

The particles that have collided with the piston no longer make a Poisson process, 
hence their distribution is much harder to control. This is our main trouble. On the 
other hand, we will be satisfied with much weaker estimates than those in the previous 
section. Also, many arguments and constructions in this section are similar to those in 
Section 3, and we omit some details. We will focus on new ideas. 

Here our analysis is always restricted to the configurations a; e Jig. Later on we will 
put additional requirements on uj. 

Recall that for uj & Qq the piston velocity is small, |V(t)| < Beq, sec (3.44), on the 
zero recoUision interval (0,Ti). Hence, the velocities of gas particles that experience one 
collision with the piston on the interval (0,Ti) are bounded 



Z'l{t,uj)\<k\Zl{t,u)\ + 



C InL 
L 




and 




^'l.min <\V\ < Vl.max 



(4.1) 
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with 

^^l,min := ^^min " 2-BSo and Vi^^ax := ^^max + 2-B£o (4.2) 

The first time of the second recolhsion T2 = 22(0;) is defined by T2 — sup^^ol^'a = 0}, 
see (3.7). Due to (3.44) and (4.1) the following bound can be easily obtained as in the 
proof of Lemma 3.2: 

rr. L L 3L , 
T2< + < (4.3) 

^max ^max ^J^^O ^max 

Now let {x,v) E and {xt,Vt) = F^{x,v) for t > 0. Denote by 

51(0;,^) = min{t : Xt — X{t)} 

the time of the first collision with the piston. The region 

Gt{t) := {{xt,yt) e Gt{t) : siix,v) < T,} 

is occupied by points that by the time t have experienced one collision with the piston, 
which occurred before time Ti. By removing the superscript + in the above formula we 
define G^{t). Let < T2 be the earliest time the piston interacts with the particles 

ix,v)e[Gtit)\Gtit)]uGtit) 

The time T^, is a random analogue of introduced in Section 2. During the interval 
(Ti,T^) the piston only interacts with the particles from Gq (t) U G^{t), hence their 
velocities must be bounded by (4.1). Denote by 

Xi{t) = {{X,V) : X = X{t) + 0, -Vl.max <V < -fl,min} 

U{{X, V) : X^ X{t) - 0, Vi^^in <V < Vi,m^} (4.4) 

two immediate one-sided vicinities of the piston which contain all "incoming" particles 
for every t e {Ti,n). 

Again, as in the previous section, we define a subinterval (Ti, 5*2) C (Ti, T*) on which 
the piston is slow enough: 

Definition Let (Ti, S2) C (Ti, T*) be the maximal time interval during which 

(a) \V{t)\ < Vi 

,min) 

(b) X,{t)GGo{t)UG,{t). 

The condition (b) means that the particles with velocities fi,min < I'^^l < i^i,max that 
are about to interact with the piston at time t have interacted with the piston during 
the interval {0,t) at most once, and if they did, the interaction occurred before Ti. 

Next we estimate how large the interval (Ti, 5*2) is. Suppose 

\V{t) -W{t/L)\ < A (4.5) 
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for all t < S2 and some small A (we will later estimate A and show that A — > as 
L — > 00). This immediately imphes \V{t)\ < A + Bsq, according to (2.26). Integrating 
(4.5) with respect to t gives 

\X{t) - LY{t/L) \ <tA 

and hence 

\X{t)-L/2\<{A + Bso)t (4.6) 

Proposition 4.1 // (4-5) holds for t < S2 with some small A > 0, then T2 — S2 < 
CL{A + Eq), where C > is a constant. 

Proof. If S2 — T2, then the statement is trivial. If 5*2 < 72, then either S2 <T^, and so 
the condition (a) or (b) in the previous definition fails at time S2, or S2 — < T2. Note 
that the condition (a) cannot fail abruptly, since we assume \ V(t)\ < A + Beq, on (0,6*2), 
i.e. V{t) remains small. If (b) fails, then at time 5*2 the piston collides with a point (,x, v) 
such that fi,min < 1"^! < 'i^i,max and the past trajectory {xt,Vt) := F^~^'^{x,v) of that point 
for t G (0, 5*2) hits the piston at some time ti > Ti. If S2 = < T2, then at time 5*2 the 
piston recoUides with a gas particle {x,v) whose past trajectory {xt,Vt) :— F^~^'^{x,v) 
experiences the first collision with the piston at some time ti > Ti. In the last case, by 
(4.5) 

\V\ < t^max + 2\V{ti)\ < i^niax + 2(A + BEo) 

In either of the above two cases, the trajectory {xt,Vt) collides with the piston twice - 
once at time ti > Ti and the second time at ^2 = 'S'2. Denote by Xi,X2 the positions 
of the piston and by Vi,V2 its velocities at times ^1,^2, respectively. Without loss of 
generality, assume that our trajectory {xt,Vt) lies to the right of the piston. Note that 
the speed \v\ — \vt\ ior ti < t < t2 satisfies 

\v\ < max{t;i,max, ^^max + 2(A + BEq)} 

< v^^ + 2{A + BiEo) (4.7) 
with Bi — max{B,B}. Then we write an obvious identity 

\v\{t2-t,) = {L-X,) + {L-X2) 

hence, by (4.6) and (4.7) 

^ _^ ^ I-(A + g,-o)(f,+t2) 
' v^^ + 2{A + BiEo) 

On the other hand, consider the particle that experiences the very first recoUision with 
the piston (this happens at time Ti). After the collision, that particle acquires velocity 
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\v{Ti + 0)1 > Vraax — '^\V{Ti)\. The next collision of this particle with the piston occurs 
after T2. Therefore, 

T,-T^ < V2 + \X{T ,) - L/2| + L/2 + iXj^) - L/2\ 



< 



w-2|y(ri) 

L+{A + Beo){T,+n) 



2(A + Beo) 
Combining the above estimates gives 

T2-S2< (T2 - Ti) - {t2 - h) < CL{A + £0) 

with some C > determined by B and B. □ 

Next, we study the dynamics during the time interval (Ti,5'2). We again define 
the density of colliding particles "on the piston" q{v,t) by the equation (3.9) and the 
functions Qi{t), i = 0,1,2, by (3.10)-(3.12). We emphasize that now, unlike what we 
had in the previous section, the density p{x,v,t) essentially depends on uj (at least for 
{x,v) G G+(t)), hence q{v,t) and Qi{t) will depend on u) not only through the piston 
position X{t) but also through the surrounding density p{x, v,t). 

In the previous sections we also introduced the deterministic density p{x,v,t). So, 
now we can define the corresponding deterministic density "on the piston" 

pix{t) + o,v,t) a v<o 
^^''^*^-\p(x{t)-o,v,t) a v>o ^^-^^ 

cf. (2.17) and the deterministic functions Qi{t), i = 0, 1,2 by the equations similar to 
(3.10)-(3.12) but with tildes over the corresponding functions. We use tildes to distin- 
guish these deterministic functions from the random ones. 

We now compare the random functions p, q, Qi with their deterministic counterparts 
on the interval (Ti,S'2). Since the transformations : '■ ~^ 

G~^(t) (recall that G~^{t) = F*(G+)) are invertible area-preserving maps, then so is the 
map 

^* = F* o F~* 

which takes G^{t) onto G'^{t). The next lemma easily follows from Theorem 3.8: 

Lemma 4.2 Let t < S2 and uj G Oq- ^^r every {x,v) e G^{t) put {x,v) :— **(a;, v). 
Then 

\x — x\<C\nL and \v — v\ < CL^^ InL 
with a constant G > 0. We also have 

p(x,v,t) ^p{x,v,t) 
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Next, the properties of the density p{x, v, t) stated in Lemma 3.3 obviously hold for 
the deterministic density p{x, v, t) on the entire region G~^. That is, the density p{x, v, t) is 
piecewise smooth with \dp{x^v^t) / dx\ < D'l/L and the discontinuity lines of 
have slope of order 0(1/L) (so they are almost parallel to the x axis). Therefore, as in 
(3.13), we have 



dX 



< ^ (4.9) 



Lemma 4.2 allows us to compare Qi and Qi considered as functions of X in the following 
way: 

Lemma 4.3 Let t < S2 and oj eVt*Q. Then for i = 0, 1, 2 

\Qi{t) - Qi{t)\ <CL-HnL 

Note that Qi{t) is defined through q{v,t) which uses the (random) position of the 
piston X{t) for the given a;, see (4.8). The above lemma shows that the dependence of 
Qi{t) on u through the density p{x, v, t) (which itself depends on ou) is very weak, because 
InL is small. In other words, the density p{x, v, t) of the gas surrounding the piston 
fluctuates with u very little. 

The following theorem is an analogue of Theorem 3.4. 

Theorem 4.4 For all sufficiently large L there is a set fll C fl^ of initial configurations 

of particles such that 

(i) there is a constant c > such that 

P{ni)>l-L-^^''^''^ (4.10) 

(ii) for each configuration cu & fl^, for each time interval 

{t,t + At) C {Ti,S2) 

such that 

the change of the velocity of the piston satisfies 

V{t + At)-V{t)^V{t)At + x (4.12) 

where 

T){t) = a[Q^{t)V''{t) - 2Qi{t)V{t) + Q2{t)] (4.13) 

and 

lxl<ciH^ (4.14) 

with some global constant C > 0. 
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Remark. Note that our bound (4.14) on random fluctuations represented by x is niuch 
weaker than (3.17) in Theorem 3.4. This is due to the lack of a good control over large 
deviations for the distribution of gas particles, as it will be clear from the proof. 

Proof. Our argument basically goes along the lines of the proof of Theorem 3.4. But it 
involves a good deal of new constructions, which we describe in detail. The first step is 
the velocity decomposition scheme, see (3.27), 

k 

V{t + At) = (1 - £k)V{t) +eJ2vj + x^'^ + X^'^ (4.15) 
The error terms x^^^ X^"^^ defined in Section 3 after (3.27), and they are bounded 

by 

for = 1, 2, see (3.32). We will see later that 

k < const • At (4.16) 

hence 

I X^"^ I < const- (At) ^ (4.17) 

for u^l,2. 

Next, we need a crude upper bound SV on possible fluctuations of the piston velocity 
V{s) during the interval {t, t + At), as defined by (3.28). One can be easily derived from 
(4.15): 

SV < 2v^^^ek < const • At (4.18) 

We now define two regions Di = Df U Di and D2 = D2 U D2 on the x, v plane 
by equations (3.29)-(3.30), where fmin and t^max must be replaced by vi,min and i'i,max, 
respectively. As it is explained in Section 3, all the particles in Di (at time t) will 
necessarily collide with the piston during the interval (t,t + At). And all the particle 
that actually collide with the piston during that interval are contained in D2. Therefore, 
we have 

kf <k'^ < kf (4.19) 

where k^,kf,kf are defined around (3.34). This gives upper and lower bounds on the 
number of colliding particles. 

Our next step is to estimate the numbers k^, i — 1,2. Since kf no longer have 
Poisson distribution, we estimate them by using a new approach. Let D be one of the 
four trapezoids Df, i = 1,2, and let kf)^^ be the (random) number of particles in D at 
time t. Obviously, kD,uj is equal to the number of particle in F^*{D) at time zero. 

We note that the trapezoid D has height t'i,max — 't'l.mm = const and width 0(At), 
hence its area is bounded by 

diAt <\D\<d2At (4.20) 
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with some constants < di < ^2 < oo- 

Let us examine the region F~*{D) more closely. Since (t, t + At) C (Ti, 5*2), it follows 
from the condition (b) in the definition of S2 that D C Go(^)UG'*(t). Put Dq = Dr\Go{t) 
and Di = Dn G^{t), then F-\D) = F-\Dq) U F-\Di). To emphasize the dependence 
of the flow F* on uj we will write F^ for F*. Now the part F^\Do) wiU be obtained 
by a simple linear transformation of Dq without collisions with the piston, hence it will 
be actually independent of 00. The part F~^{Di) is also obtained by pulling the domain 
Di back in time, but one collision with the piston will occur along the way. Since the 
position and velocity of the piston are random, then the domain F~^[Di) will depend on 
u. Hence, the domain := F^\D) wiU depend on uj. The initial number of particles, 
kD,ui, in a randomly selected domain D,^ certainly need not be a Poisson random variable. 

To estimate ku^uj we fix some At satisfying (4.11) and construct finitely many domains 
Dn C G, 1 < n < A^*, which have the following property. For every G fig, every 
t e (Ti, 5*2 — At), and every trapezoid D defined above, there are two domains D^" 
such that 

D^, C FJ*(L>) C (4.21) 
and the area of the difference is relatively small: 

Pn"\^n'| <X^'^Pu,| (4.22) 

with some x^^^ — as L — > 00. We say that D^' and D^" approximate "from inside" 
and "from outside", respectively. We denote the collection of the domains D„, for the 
given At, by 

C = CAt = {Dn}n=l 

We postpone the construction of D^s for the moment and derive immediate benefits of 
the above approximation. The inclusion (4.21) implies 

k„',cj < kD,uj < kn",ij (4.23) 

where kn,ui is the (random) number of particles in the domain at time zero. Since 
Dn, for each n, is fixed (independent of cu), the random variable kn,ui does have a Poisson 
distribution with mean value 

jDn 

According to Lemma A. 2, we have 

P (a; : \kn,u; - Xn\ > Bsl)Oj < 2e-^^' (4.24) 

for any B < b\Ov^, where 6 > is a constant and c = c{b) > another constant. We will 
specify the value of B — BAt (one for all DnS in Ca*) later. Due to (4.20), it is enough 
to require 

BAt < (4.25) 
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Then we define ^^(Ai) as the set of configurations a; e Jig satisfying 

\kn,c. - K\ < Bm\Ik for 1 < n < iV* (4.26) 
Then by (4.24) we have 

P{% \ ni{At)) < 2A^,e-^-^i* (4.27) 
Next, for all u; e flKAt) the bounds (4.23) and (4.26) imply 



An' — BAty < kcul < A„" + BAty ^n" (4.28) 
Furthermore, consider the quantity 

^Du> — L^ / p(x,v,Q) dx dv = / pix,vA)dxdv 
Jf-\d) Jd 

Due to the inclusion (4.21) we have 

An' < )^D,ijj < An" 

It follows from (4.22) that 

(1 - cx('^)Ao,, < An' < An" < (1 + cx('))A^,, (4.29) 

for some constant c > (determined by TTmax, TTmin, vi, and V2 in (P3)). 
Consider the deterministic quantity 

\n = L^I p(x,v,t)dxdv 
Jd 

It easily follows from Lemma 4.2 and the properties of the function p(x, t) that 

\\D,^-~\D\<CL-^\nL\D (4.30) 
with some constant C > 0. Combining (4.29) and (4.30) we arrive at 

(1 - X^'))Ad < An' < A„" < (1 + x('^)Ai5 (4.31) 

with 



X 



(4) 



const ■ (x^^^ +i^"MnL) (4.32) 

The bounds (4.28) and (4.31) will give the desired estimate on the number k^^ui- 

We now construct the domains that approximate the domains D^^ = F^^{D) for 

all a; e Qq trapezoids D defined above. We first fix a; e and a trapezoid D and 

will construct two special domains D' , D" that approximate Di_^ — F~^(D) from inside 

and from outside, i.e. such that D' C C D". 

The domain is obtained by pulling D back in time. We consider its "trajectory" 

D~ = F^~^{D) for < s < t, so that Dq = and = D. We examine the shape of 
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the domain D~ and how it changes as s runs from t down to 0. Recall that the trapezoid 
D is adjacent to the piston at time t (it is about to run into the piston at that time). 
As s goes from t downward, the domain D~ comes off the piston and travels to a wall 
(as in a movie running backward). During that period, the map F~^^~'^^ restricted to D 
is linear, hence D~ is still a trapezoid. But since the velocities of points {x, f ) G -D vary 
(from fi^min to fi^max), the trapezoid will be "skewed" - its "outer" edge |f | = fi,max 
will move toward the wall faster than the other edge |f | = fi,min- By the time it reaches 
the wall, D~ will be a long slanted trapezoid stretched the distance 0{L) along the x 
axis. Every vertical line (parallel to the v axis) will intersect D~ in a segment of length 
0{At/L). As s runs farther down, a collision with the wall occurs, and a new part of D~ 
appears directly across the x axis, moving now toward the piston. Its shape will be also 
that of a long narrow trapezoid, whose vertical "thickness" is 0{At/L). Eventually it will 
move all the way (the distance 0{L)) from the wall to the piston and contact the piston 
at some time s^: = s^{D), see Fig. 5. Note that so far is completely independent of 




Figure 5: The domain . The arrows show its motion as s decreases. 
The "outgoing" part of D~ , as it comes off the piston, is shown on the right. 



After the collision with the piston, a new part of appears across the x axis, coming 
off the piston and moving back to the wall. This is the most interesting part, it will be 
determined by the piston position and velocity, hence it will actually depend on uj. Recall 
that the part of D~ running into the piston before the collision (the "incoming" part) 
is a narrow trapezoid bounded by two almost horizontal lines with slope 0{1/L) lying 
a distance 0{At/L) apart. We call them the upper and lower boundaries of Dj . The 
part of coming off the piston after the collision (the "outgoing" part) also has upper 
and lower boundaries (due to the reversal of velocities at collision, though, the upper 
boundary after the collision is the image of the lower boundary before the collision, and 
vice versa), see Fig. 5. The boundaries of the outgoing part of D~ are quite irregular 



55 



and depend on u. It is them who need careful approximation. 

The shape of the outgoing boundary of D~ will be determined by the piston velocity 
V{u) during the interval s < -u < s*. Since the collision between D~ and the piston occurs 
during the zero-recoUision interval (s* < Ti), we can apply the results of Section 3. By 
Theorem 3.7 the piston velocity behaves as a Holder continuous^ function: 

\V{u + h)- V{u) I < CL-^ In LVh^: dh (4.33) 

for all h e (L~^,l). We will only consider h satisfying L~"i < h < L""^ -^th some 
< a2 < ai < 2. Note that in this case 

Ldh/h oo and dh/h (4.34) 

as L ^ oo. 

Let us consider the upper boundary of the outgoing domain D~ after the collision. 
It is a continuous curve that can be parameterized by the collision time (or contact 
time between Dj and the piston), call it u, and then it becomes a parametric curve 
{xs{u),Vs{u)) with parameter u e (s,s*). This means, in particular, that 

F'^{Xs{u),Vsiu)) = {Xs+w{u),Vs+wiu)) 

ior s < u < and w <u — s, and the point 

F''-'{xs{u),Vs{u)) = {xu{u),Vu{u)) 

is the endpoint of the upper boundary of D~, i.e. Xu{u) — X{u) is the piston coordinate 
at time u. 

Claim 1. For each u e (0, s* — h) and all sufficiently large L we have 

Ixu+hiu + h) - Xu{u)\ < Beoh, (4.35) 

and 

\vu+h{u + h) - Vu{u)\ < 3dh (4.36) 

Proof. The first bound means that \X{u + h) — X[u)\ < Be^h and follows from (3.44). 
To prove the second, we use the collision rule 

Vu{u) = -v-{u) + 2V{u) 

where v~{u) is the v coordinate of the lower boundary of the "incoming" part of Du 
where it contacts the piston (at time u). Similarly, 

v^+h{u + h)^ -v-+hiu + h) + 2V{u + h) 

■'For any (i > 0, d-dense sets in the space of Holder continuous functions were constructed by Kol- 
mogorov and Tihomirov [KT], our constructions here are in the same spirit. 
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Since the lower boundary of the "incoming" part of is a straight hne with slope 
0(1/L), and the piston velocity is bounded, then 

\v~j^h{u + h) — v~{u)\ < const ■ h/L < dh 

where the last bound follows from(4.34). This and (4.33) imply (4.36). □ 
We now fix some h satisfying 

{At/iy/^ <h< 2(At/L)i/2 ^437) 

With our restrictions (4.11) on At, this imphes L~^/^ < h < L~^/^, hence (4.34) will 
hold. Also, we note that 

dh/{At/L) ^0 as L ^ oo (4.38) 

hence dh will be much smaller than the "thickness" of the domain D~ . 

Next we put a lattice on the x, v plane with the x-spacing rh and the v-spacing rdh, 
where r > is a sufficiently small constant (for example, r = S£o/10). The lattice sites 
are 

{xi,Vi) ^ {rhi,rdhj), e ^ (4.39) 

Now we are ready to define the boundary of D' C Dq = F~^{D) approximating 
Dq from inside. We start with the upper boundary of the "outgoing" part of Dq . It 
is parameterized as {xo{u),Vo{u)), see above, with < -u < s*. We consider discrete 
parameter values = kh, k = 0,1, . . . , [s^/h]. For each such k, consider the point 

F''''{xo{kh),vo{kh)) = {xkh{kh),Vkh{kh)) 

This is the endpoint of the upper boundary of D^f^, so that Xkhikh) = X{kh) is the 
piston coordinate at time kh. Now we pick a site of the lattice (4.39) closest to the point 

{xkh{kh),Vkh{kh)) and lying on the same side of the piston as the domain D'^y^. Call this 
site {xl,yl). Next, we adjust this site by moving it down (along the v axis) the fixed 
distance lOdh and obtain the adjusted point (a;^,^^): 

xt = xl, vt ^vl- 104 (4.40) 

We note that Claim 1 (along with the smallness of r) implies 

\4 - = 14 - 4+il < '^Bsoh (4.41) 

and 

K - vt+^\ = K - vl+,\ < Mn (4.42) 

Next, let 

ixlvl) := F-'\xt,vt) 
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Now join the points {x1,v^) and t'^+i) with a straight hne segment, call it L^. 

Then the upper boundary of the domain D' is made by the segments L^, i.e. it is 

Ufc=o 

Claim 2. The upper boundary of D' lies completely inside the domain Dq — F~^{D). 

Proof. It is enough to show that each link lies inside Dq . It is easy to see that the 
map F^^'' restricted to the link Lk is a linear map, and the image :— F^^[Lk) is a 
straight line segment joining the points (a;^,^^) and 

We denote by 

the lengths of the projections of onto the x and v axes, respectively. Note that 
Vi,mm — Beq < It'fe+il < ■i'i,max + Beq. Hcncc (4.35) implies that 

(^^i,min - 2Beo)h < \Ll\^ < (vi,niax + 2Beo)h 

i.e. is of order h. By (4.36), |L^|^ < Sdh- Hence, the slope of LI is 0{dh/h) ^ by 

(4.34). Therefore, is almost a horizontal segment. The inequality (4.36) shows that 
the upper boundary {xkh{u),Vkh{u)) of for kh < u < kh + h, which lies directly above 
L^, oscillates in the v direction by less than 3dh. The adjustment (4.40) then ensures 
that the segment lies entirely below the upper boundary of the domain D^^. Also, 
(4.38) implies that the domain D^f^ is much "thicker" than the distance between its upper 
boundary and L^, hence lies entirely inside that domain. Hence, — F~^'^{L%) lies 
entirely inside Dq = F~'''^{D'i^f^), proving our claim. □ 

It is also clear from the above argument that the segment lies the distance < 20dh 
below the actual upper boundary of the domain Z^^^^ so that the area between them is 
bounded by const- hdh. Hence, the total area between the upper boundaries of Dq and 
D' is bounded by 

const • hdh ■ (s*/ h) < const • Ldh 

In the same way we construct the lower boundary of the new domain D'. The only 
difference is that we adjust the selected sites of the lattice (4.39) by moving them up, 
so that the joining segments will be again inside Dq . The upper and lower parts of the 
boundary will give us a new domain approximating the "outgoing" part of Dq after the 
collision with the piston, i.e. the part F^^{D fl G^{t)). Note that this part may itself 
experience one collision with the wall (during the interval (0, s^)), then it will consist of 
two connected components lying across the x axis. 

It remains to approximate the "good" part of F^^{D) = Dq , which has not interacted 
with the piston, i.e. the part F^^{D fl G'o(t)). That one, as described above, consists of 
one or two trapezoids bounded by two almost horizontal straight lines a distance 0{At/ L) 
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apart. This is an easy task. We simply replace the upper (lower) boundary of the domain 
F~*{D n Go(^)) with a polygonal line joining some sites of the lattice (4.39), picked one 
on each vertical line x = Xi crossing the domain D~ and lying inside (outside) of this 
domain and no farther than dh from the original boundary (recall that r in (4.39) is very 
small, so the above choice is possible). That gives the boundary of D' . Note that if we 
number the selected sites consecutively (say, from left to right), and call them (x|,y|), 
then the neighboring sites will again satisfy (4.41)-(4.42). 

The complicated construction described above produces a domain D' C D,^ — F~^{D) 
such that 

l^c. \ ^1 < const • dhL = const -InlVh (4.43) 

In a completely similar way we construct another domain D" that approximates D^^, from 
outside and satisfies 

\D"\D^\< const • dhL = const ■ In L^/h (4.44) 

(we just need to adjust the sites along the upper boundary by moving them up and the 
sites along the lower boundary by moving them down, so that the new boundaries will 
be completely outside of D^^). Since \D\ — 0{At), the inequalities (4.43)-(4.44) imply 

\D"\D'\< const- In LVh{At)-^\D\ 

hence we get (4.22) with 

X^^^ = const ■ In iVh (At)"^ (4.45) 

Note that x^^"* ^ as L — > oo due to our choice of h in (4.37) and restrictions on At in 
(4.11). 

Our domains D', D" are constructed around D^^, hence they depend on u. However, 
the boundaries of D', D" are polygonal hues whose vertices are derived from the sites of 
a fixed lattice (4.39). Therefore, the total number of distinct domains D',D" is finite. 
Now we estimate their number A^*. 

First, the total number of the sites of the lattice (4.39) in the relevant area < x < L, 

\v\ < fl,max is 

Ki = const • (L/h) ■ {2v^^/dh) < const • L {hdh)~^ (4.46) 

Recall that the domain Dq — F~*{D) consists of at most four connected components. 
The part F~*{D fl G'*(t)), or the "outgoing" part of Dq created after the collision with 
the piston, consists of at most two components. And the part F~^{D r\Go{t)) consists of 
at most two components, each is a regular, trapezoidal region. The upper (lower) part 
of the boundary of D' {D") will then consists of at most four disjoint polygonal lines 
constructed from the lattice sites (4.39). In each polygonal line, the first point can be 
constructed from any of the Ki sites of the lattice, see (4.46). But other points can be 
constructed consecutively, and each point must be constructed from a site selected from 

Kq := {ABeo/r) ■ (8/r) = const (4.47) 
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nearest neighbors of the previously selected site, according to (4.41)-(4.42). It is also 
clear that each polygonal hne has at most 

M := const • L/h 

links (vertices). Therefore, the total number of ways to construct one polygonal line does 
not exceed KiK^^ . The total number of ways to construct the entire boundary of D' 
{D") is then less than KfR^^. This gives an upper bound on the number of distinct 
domains D' {D"): 

N, < KlKl^^ < const • L^{hdh)-^ ■ e^°°^*-^/'^ (4.48) 

Of course, the exponential factor is dominant here and can absorb all the others. 
We now set the value of BAt in (4.25)-(4.28): 

BAt = C7 \/L/h (4.49) 
with a sufficiently small constant cy > 0. Then (4.27) and (4.48) imply 

P{n*o \ fll{At)) < const • e-^^^/^ (4.50) 
with some constant Cg > 0. Recall that h satisfies (4.37), hence 

P(Q* \ ni{At)) < const • e-^°°«t-^'^'(^*)"'^' 
Next, since At satisfies (4.11), we have 

P(Q* \ ni{At)) < const • e-^"-^^*--^"^' (4.51) 

Also, by (4.49) and (4.37) 

BAt < c,L^'\At)-''' 

which implies (4.25) since At ^ L~^/^. 

Next, for each u e Q.\{At), each i e (Ti, S2 — At), and any trapezoid D defined above, 
the number of particles kD,uj in D satisfies (4.28)-(4.32) with x^^^ given by (4.45) and h 
given by (4.37), hence 

kD,u} = Ad + XD,a; 

with 

\XdA < const- (^B At L\/At + In LVh + LlnL^ 

< const- L"HnL{Atf'^ 

We note that = O(L^At), hence \xd,u)\ A^ for all At satisfying (4.11). By the 
way, this easily implies a rough bound kD,ui < 2Ad <const-L^ At, which justifies (4.16) 
and hence (4.18). 
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Next, using (4.19) and applying the above estimates to each oi kf, i — 1,2, gives 

k^L^l p{x,v,t)dxdv + x' 

JDi 



with 



\X 



'\ < const • [L^/^ In L (Ai)^/^ + 5V At] 



In a similar way we can estimate the other random factor in the main decomposition 
formula (4.15), which is 

k 

Z = E^i 

and get 

Z^L^f vPix,v,t)dxdv + x" 

JDi 

with the same upper bound on x" thai on x' above. 

Using the smoothness properties of the function p(a;, v, t) described around (4.9) gives, 
cf. (3.38)-(3.40), 

/ p{x, V, t) dx dv = {Qi{t) - Qo{t)) At + x^^^ 

JDi 

[ vp{x,v,t)dxdv= {Q2{t) - Qi{t)) At + x^^^ 
|X^"^| < const • \{Atf/L + 6VAt 



and 



with 



for u = 4,5. 

Therefore, we get 

V{t + At) - V{t) = Vit) At + x (4.52) 

where 

Ixl < const - [L-^/^\nL{Aty/'' + 5VAt+{Atf] (4.53) 

Recall that 5V = 0{At) by (4.18). It is now easy to check that the second and the 
third terms are much smaller than the first one. This proves (4.12)-(4.14). It does not 
prove Theorem 4.4 yet, because we have fixed one (arbitrary) value of At, and our set 

Ql depended on At. 

Actually, for our main purpose it is enough to prove Theorem 4.4 for just one value of 
At, namely for At = 1/L^/^, as we will see later. But at a little extra effort we can prove 
our theorem for all At satisfying (4.11). We do that next. Divide the interval (4.11) 
into subintervals of length e~^. That is, fix a finite collection of points (Ai)„ = ne~^ for 
n — rii, . . . ,n2 with rii — L~^/^(ln L)^e^ and n2 — L~^l'^e^. Then we define 
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The bound (4.51) then imphes 



p{n*o \ ni) < const • 6^-^°"^^*-^"^' 

which is obviously sufficient to maintain the bound (4.10). 

Now for any At satisfying (4.11) we find (At)„ such that \At — (At)„| < e~^. For 
any uj E G r2^((At)„) we have all the above estimates with At replaced by (At)„. 
This replacement only causes an exponentially small error, e~^, in our estimates. It 
will not spoil our bounds, which are all polynomial in L. This completes the proof of 
Theorem 4.4. □ 

The next theorem is an analogue of Theorem 3.5. 

Theorem 4.5 Assume that Eq > is small enough. For all sufficiently large L, for each 
configuration u) eVL\ and for all t e (Ti, ^2) we have 

(i) there is a constant B > such that 

\V{t)\<Beo (4.54) 

(ii) there is a constant Co > such that 

\V{t) - Vo{t)\ < CoL-^/'^lnL (4.55) 

where Vo{t) is defined by 



Qo{t) 



whenever Qo{t) 7^ and by 



otherwise. 



Proof of this theorem very much repeats that of Theorem 3.5. The first half of it, up to 
the formula (3.59) can be copied almost verbatim, with only replacement of Si by S2, Qi 
by Qi, and V by V. We omit that part. The rest of the proof requires more substantial 
modifications, and we give it in detail. 

First, recall that Theorem 3.5 deals with t G (0,6*1). On the interval (0,6*1) the 
functions Qi{t), i = 0, 1,2, are independent of uj, they are defined by equations (3.9)- 
(3.12) where p{x,v,t) was in fact the deterministic density now denoted by p{x,v,t). 
Therefore, our functions Qi{t) for t > Ti, are natural continuations of Qi(t) beyond the 
interval (0,6*1) = (OjTi) and they have the same properties, cf. (3.13) and (4.9). Hence, 
the function Vo{t) defined by (4.56)-(4.57) for t > Ti is the continuation of Vo{t) defined 
by (3.46)-(3.47) on the interval (0,Ti). 
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Next, in the first half of the proof (which we omitted since it almost coincided with 
that of Theorem 3.5), we must introduce < S2 as the first time when (4.54) fails. Now 
we prove (4.55) for all t <U with some constant Co > (independent of the choice of B 
in (4.54), which is to be made yet). 

The bound (3.45) proved for all t < Ti imphes that (4.55) holds for at least some 
t > Ti. Next, if (4.55) fails for any t < t^, then let t e {Ti,t^) be the first time (4.55) 
fails. Denote by 

the maximal allowed time increment in Theorem 4.4. Let s — t — Ag. Due to Theorem 4.4 

V{t)^Vis)+V{s)Ao + x (4.58) 

with 



< CL-i/^ In L (Ao)^/^' = In L 

Due to the analogue of (3.58) obtained in the first half of the proof of the theorem, 

Vo{t)^Vo{s)+Xo (4.59) 

with 



< 



Eq Eq Aq Eq Eq 



L 

For brevity, put U{s) — V{s) — Vo(s) for all s. Subtracting (4.59) from (4.58) then gives 

Uit)^U{s)+V{s)Ao + x' (4.60) 

with x' = X — Xo; so that for large L 

Ix'l < 2CL-2/^lnL (4.61) 

Now assume, without loss of generality, that U{t) > 0. Since (4.55) fails at time t, we 
have 

U{t) >CoL-^/^ In L (4.62) 
Now consider two cases. If U{s) < 0, then by the analogue of (3.59) 

U{t) < \V{s)\Ao + 1x1 < E2 \U{s)\Ao + Ix'l < L-i/^lnL 

for large L, which contradicts to (4.62). If U{s) > 0, then, again due to (4.60) and the 
analogue of (3.59), 

Uit)<Uis)[l-E,Ao]+x', 



hence 



U{s) > ^^^^>{U{t)-x'){l + E^Ao) 

> U{t) + U{t) El Ao-2x' (4.63) 
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Now, if Co in (4.55) is large enough, say Co = ?>C/Ei, then U{t)EiAo > 2x' by (4.62) 
and (4.61). This fact and (4.63) imply U{s) > U{t), so (4.55) fails at an earlier time 
s < t, a contradiction. Hence, (4.55) is proved for all t <t^ and Cq = 3C/Ei. 

Lastly, the remaining part of the proof of Theorem 3.5 can be repeated verbatim, 
concluding the proof of Theorem 4.5. □ 

We finally prove the convergence, as L — > cxo, of the random trajectory of the piston 
to the solution Y{t/L), W{t/L) of the hydrodynamical equations described in Section 2 
for all t e {Ti,S2). 

Theorem 4.6 Assume that Sq > in (P5) is small enough. Then, for all large L and 
all cu & 0,1, there is a constant C > such that 

|y^(r,a;)-y(r)|<^ (4.64) 

and 

\W,{T,u;)-W{r)\<^^ (4.65) 
for all min{Ti, Ti/L} < r < min{T2, S2/L}. We also have 

T2-S2< CeoL (4.66) 

and 

T2 - S2/L < Ceo (4.67) 
The fluctuations of the function S2 — S2 {uj) for u & fll are bounded by 

sup \S2{u;)-S2{u')\<^^ (4.68) 



Proof. According to (2.29), the deterministic function Y{t) satisfies 

dY{T)/dT - F(Y",r), y(0) = 1/2 (4.69) 
Now Theorems 3.7 and 4.5 imply that for all u the random trajectory satisfies 

dYL{T,uj)/dT = F(r,r) + x(t,^), Yl{Q,uj) = 1/2 (4.70) 

with some C \ L 

|x(T,a;)| < -jTp- 

Recallthat \dF{Y,T)/dY\ < k, see (2.30). Therefore, the difference ^^(r, cu) := Yl{t,u)- 
Y{t) satisfies 

I / \ I i ^ / \ I G 111 ^ 

|Zi(r,a;)|<«:|Zi(r,a;)| + -^ 
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and Zl{0, cu) — 0. Using the standard Gronwall inequality in differential equations, see, 
e.g.. Lemma 2.1 in [TVS], gives 

and 

\Z',ir,u)\<^e- 

for all r < S2/L, which imply (4.64) and (4.65). 

Next, (4.65) enables us to apply Proposition 4.1 with A = CL~^/^lnL and thus prove 

(4.66) . Now we employ the same argument as in the proof of Theorem 3.8. By (4.65), 
random fluctuations of the piston velocity are bounded by CL~^/^ In L. Hence, random 
fluctuations of the velocities of particles that have had one or two collisions with the 
piston are bounded by 4CL~^/^lnL. The random fluctuations of the positions of both 
piston and particles at every moment of time t < min{T2L, S2} are bounded by the same 
quantities (with, possibly, a different value of C) in the coordinate y = x/L. This implies 
(4.68) in the same way, as (3.92) in Theorem 3.8. Now we can apply Proposition 4.1 to 
the deterministic dynamics constructed in Section 2 with the same result, and thus prove 

(4.67) . □ 

The bound (4.68) shows that the limit 

T** := lim S'2(cu)/L 

does not depend on G il*. We have |t** — 2/vjnax.\ < const • Eq, due to (4.67) and the 
estimates in Lemma 2.9. 

The convergence claimed in Theorem 1.2 is now proved on the interval (0,r**), but, 
generally, r** < with — r** = 0{6o), so we may still be 0{so) short of our target value 
T*. To extend our results all the way to r* we need to redeflne S2 and the neighborhood 
Xi introduced by (4.4) more accurately. We need to set 

A'i(t) = {{x, v) G G+{t) : X = X{t) + 0, v < 0} 

U{(x, v) eG+{t): x = X(t) - 0, v>Q} (4.71) 

The domain G^{t) = F^[G~^) is random (it depends of u), i.e. we should write G^{t) = 
G^{t,uj), and this is why we could not adopt the above definition of Xi{t) earlier and 
opted for a cruder one (4.4). But now, as we have just shown in the proof of Theorem 4.6, 
random fluctuations of the velocities and positions (in the y coordinate) of the particles 
and the piston are bounded by const • L~^/^lnL. Hence, at every time moment t < 
min{r2-C/, 5*2} all the domains G~^{t,uj), are close to each other - the distance between 
G~^{t,uj) and G~^{t,uj') for ui,ui' G in the Hausdorff metric on the coordinate plane 
y,v is bounded by const • L~^/^lnL. By the same reason, every domain G+(t,u;) will be 
0(L~^/^lnL)-close in the Hausdorff metric to the deterministic domain Q^[t/ L) deflned 
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in Section 2. Hence we can easily make G'^{t) in (4.71) independent of uj by, say, taking 
the union 

G+{t)^U^eniG+{t,u;) (4.72) 

It is important to note that all the gas particles colliding with the piston at time t for 
every G il* will be in Xi defined by (4.71)-(4.72). With this new definition of Xi 
replacing (4.4) and with 5*2 changing accordingly, we have 

lim \S2/L-n\ = 

which follows from the 0(L~^/^ In L)-closeness of the random dynamics to the determin- 
istic dynamics on the y,v plane. Thus we extend all our results to the interval (0,t*). 

5 Beyond the second recollision 

The main goal of our analysis in Sections 3 and 4 is to prove that under suitable initial 
conditions random fiuctuations in the motion of a massive piston in a closed container 
filled with an ideal gas are small and vanish in the thermodynamic limit. We are, how- 
ever, able to control those fiuctuations effectively only as long as the surrounding gas of 
particles can be described by a Poisson process, i.e. during the zero-recoUision interval 
< T < Ti. In that case the random fluctuations are bounded by const- In L, see The- 
orem 3.8. Up to the logarithmic factor, this bound is optimal, according to Theorem 1.1 
by Holley. 

During the one-recollision interval Ti < r < r^,, the situation is different. The proba- 
bility distribution of gas particles that have experienced one collision with the piston is 
no longer a Poisson process, it has intricate correlations. We are only able to show that 
random fluctuations remain bounded by L~^/^, see Theorem 4.6. Perhaps, our bound 
is far from optimal, but our numerical experiments reported below demonstrate that 
random fluctuations indeed grow during the one-recollision interval and beyond. 

At present, we do not know if our methods or results can be extended beyond the 
critical time r*, this remains an open question. We emphasize, however, that our Theo- 
rem 4.6 is the flrst rigorous treatment of the evolution of a piston in an ideal gas where 
most or all of the particles experience more than one collision with the piston. 

In order to understand what is going on beyond the critical time r*, and in particu- 
lar whether random fluctuations grow or remain small, we undertook experimental and 
heuristic studies of the piston dynamics on a large time scale. Below we describe our 
findings and discuss further research in this direction. A detailed account of our work 
can be found in [CL]. 

We set the initial density of the gas to 

My,v) = M\v\) = \^Q elsewhere ^^'^^ 
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It satisfies our requirements (P1)-(P5), in particular Vmin — 0.5 and I'max = 1, and most 
importantly Eq = 0. Therefore, by Corollary 2.13, the solution of the hydro dynamical 
equations is trivial: Y{t) = 0.5, W{r) = 0, and 7i{y,v,r) = 7ro(y,f) for all r > 0. 

To generate an initial configuration of particles, we used a random number generator 
described in [MN]. For our density (5.1), the x and v coordinates of all the particles 
are independent random variables uniformly distributed in their ranges < x < L 
and Vinin < |f I < Vraax- Our computcr program first selects the number of particles N 
according to the Poisson law with mean L^, and then generates all {xi,Vi), I < i < N, 
independently according to their uniform distributions. The parameter L changed in our 
simulations from L = 30 to L = 300. For L = 300 the system contains ^ ^ 27, 000, 000 
particles. 

Once the initial data is generated randomly, the program computes the dynamics by 
using the elastic collision rules (1.2)-(1.3). All calculations were performed in double 
precision, with coordinates and velocities of all particles stored and computed individu- 
ally. 

Figure 6 presents a typical trajectory of the piston. Here L — 100. The position and 
time are measured in hydrodynamic variables Y — X/L, < y < 1, and t — t/L. 




5 10 15 20 25 30 



Figure 6: The piston coordinate y as a function of time r. Here L — 100, N — 1000229. 

Initially, the piston barely moves about its stationary point y = 0.5. Then, at times 
T between 3 and 5, the random vibrations of the piston grow and become quite visible on 
the y-scale, but for a short while they look random, as a trajectory of a Brownian motion. 
After that the piston starts travelling back and forth along the y axis in a more regular 
manner, making excursions farther and farther away from the stationary point y = 0.5. 
Very soon, at r = Tmax ~ 8, the swinging motion of the piston reaches its maximum, 
(AF)inax = max \Yi{t) — 0.5| ~ 0.1. Then the oscillations of the piston dampen in size 
and seem to stabihze at an amphtude A ^ 0.04. At the same time the trajectory of 
the piston smoothes out and enters an oscillatory mode with a period Tper ~ 1.63. The 
velocity of the piston W{t) follows similar patterns, see Figure 7. 

Both the coordinate and velocity of the piston continue almost perfect harmonic 
oscillations for a long time with the same period Tper — 1.63 (this is independent of L) 
but the amplitudes of both Y{t) and W{t) are slowly decreasing, see Figure 8. 
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Figure 7: The piston velocity as a function of time r. The same run as in Fig. 6. 
Y 



0.54- 
0.52- 

0.5 
0.48 
0.46 






30 



35 100 



105 250 



255 900 



905 



Figure 8: The piston coordinate Y during the intervals (30,35), (100,105), (250,255), 
and (900, 905). The same run as in Fig. 6 and 7. 



The oscillations of the piston with decaying amplitude can be described, in the interval 
20 < T < 1000, approximately by 



yfr) ~ ^e-^(^-2°) si 



smo; T 



a) 



(5.2) 



with A — 0.046 and some constant A > 0. Correspondingly, W — dY/dr in the same 
interval 20 < r < 1000 is 



W{t) 



-XY{t) + Ae-^^^-^'^^uj cos uj{t - a) 
j[Q-Kr-m^_X sin uj{r - a)+u;cosu;(T - a)] 



= Aie 



-A(t-20) 



smo; T 



(5.3) 



with Ai = A\J ijJ^ + and some (5 related to a. 

To check how well our prediction (5.2) agrees with the experimental data, we com- 
puted the amplitude ^(r) as a function of time r, by fitting a sine function Y{t) = 
Asintc;(r — a) "locally", on the interval (r — 5,t + 5) for each r. Fig. 9 shows A{t) on 
the logarithmic scale, which looks almost linear on the interval 30 < r < 800. 

We used the least squares fit to estimate A = 0.00264 for the run shown on Figs. 6-9. 
Since A is small, the oscillations indeed die out very slowly. The "half-life" time (the 



68 




Figure 9: The amplitude A{t) on the logarithmic scale: experimental curve (bold) and 
a linear fit (thin). The same run as the one shown in Fig. 6, 7 and 8. 



time it takes to reduce the amplitude by a factor of two) is T1/2 = A~^ln2 ^ 263. The 
parameter A and hence T1/2 depend on the system size L. We estimated numerically that 
Ti/2 ~ L^'^, hence A ~ L~^-^. 

The key characteristics of the piston trajectory described above, in particTilar. (Ay)„iax 
W^max! ^) Tpcri appear to be independent of L. Even for L = 300 (the largest system tested 
experimentally) the piston experiences large oscillations very similar to the ones shown 
on Fig. 6 and 7. Some other quantities, such as T1/2 and the related A, depend in a 
systematic way on L. 

But most importantly, the time of the largest oscillations Tmax and the related time 
of the onset of the instability Tc, sec below, seem to slowly grow with L, very likely as 
logL. To understand this fact, we looked into the mechanism of the build-up of random 
fluctuations of the piston position and velocity displayed on Figures 6 and 7. To this end 
we plotted the histogram of the (empirical) density of gas particles in the y, v plane at 
various times < r < 30, see samples in Figure 10. The initial density (at time zero) 
is almost uniform over the domain < x < L and fmin < |f | < fmax (variations in the 
initial configuration always exist, because it is generated randomly). Then, for < r < 1, 
the piston experiences random collisions with particles and acquires a speed of order 
M~^/^ = 0(1/L), see Theorem 1.1. These small fluctuations of the piston velocity result 
in bigger changes of the velocities of the particles which leave the piston after collisions 
due to the rule (1.3). In particular, the outgoing particles on the right hand side of the 
piston have velocities in the interval (t'min + 2iy(T), t'max + 2W^(r)) while those on the left 
hand side of the piston have velocities in the interval (— I'min + 2iy(r), — Wmax + 2W{t)). 
Hence, the region in the y, v plane where the density of the particles is positive is no 
longer a rectangle with straight sides, now its boundaries are curves whose shape nearly 
repeats the graph of the randomly evolving piston velocity W{t). While the variations 
of 0(1/ L) of these boundary curves may seem small, it is crucial that on opposite sides 
of the piston they go in opposite directions. Indeed, when W{t) > 0, then the outgoing 
particles on the right hand side accelerate and those on the left hand side slow down. 
When W{t) < the opposite happens. 
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Figure 10: Six snapshots of the empirical gas density (in the x,v plane) at times r = 0, 
2.3, 4.2, 5.9, 7.4 and 18.6. 

Next, the particles that have collided with the piston travel to the wall and come back 
to the piston. Now their densities are less regular than they were initially - the regions in 
the x, V plane where the density is positive, are curvilinear domains. When they hit the 
piston, they shake it back and forth more forcefully than before, because the velocities of 
the incoming particles on the opposite sides of the piston are now negatively correlated. 
When particles on the right hand side are fast, those on the left hand side are slow, and 
vice versa. The fluctuations of the gas densities thus "cooperate" to push the piston 
harder, with a "double" force. This produces a resonance-type effect destabilizing the 
piston dramatically and the velocity of the piston W{t) experiences larger fluctuations 
than before. The velocities of the newly outgoing particles will again go up and down in 
opposite direction, on a greater scale than before. 

As time goes on, the above phenomenon repeats over and over, with larger and 
larger fluctuations of the gas and piston velocities, until the distribution of gas particles 
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completely breaks down. For L — 100, at times r ~ 10, two large clusters of particles 
are formed, one on each side of the piston. When one cluster bombards the piston, the 
other moves away from it and hits the wall, then they exchange their roles. The clusters 
have sizes of about 0.3-0.5 in the y direction and the particle velocities range from about 
0.2 to just over 1. The average velocity is about 0.5-0.6 and so the clusters hammer the 
piston periodically with period 1.6-2.0, which is close to the experimentally determined 
period of piston oscillations, sec above. 

Fig. 10 shows six snapshots of the empirical density of gas particles taken at different 
times. At r = the gas fills (almost uniformly) two rectangles {(y, v) : 0.5 < < 1, < 
y < 1}. AtT = 2.3 one can see some ripples on the boundaries of these rectangles. At 
time T = 4.2 the irregularities grow and at r = 5.9 the rectangular shape is broken down. 
Two large clusters of particles are formed, both appear in the upper half-plane f > 0, i.e. 
at that time both clusters move to the right (one toward the piston, the other away from 
it). Later the density undergoes strange formations (r = 7.4) but eventually smoothes 
out and enters a slow process of convergence to Maxwellian (r = 18.6) described below. 

The above analysis suggests that the fluctuations of the piston velocity roughly in- 
crease by a constant factor during each time interval of length one. Indeed, initial random 
ffuctuations Wa ~ 0{1/L) result in additional changes of velocities of outgoing particles 
by 2Wa- When those particles come back to the piston (in time At ~ 1), they kick its 
velocity to the level of 2Wa. Then the newly outgoing particles acquire an additional 
velocity 4Wa, etc. Over each time interval of length one the fluctuations double in size. 
This is an obvious oversimplification of the real dynamics, but it leads to a reasonable 
conjecture 

Wair) ^ ^ (5.4) 

where Wa{T) are typical fluctuations of the piston velocity at time r and C, i? > are 
constants. 
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Figure 11: The value Tc as a function of ln(l^c-^): experimental points and a hnear fit. 

We tested the above formula numerically as follows. Let Wc > be some preset 
critical value of the piston speed and = inf{r > : |Vr(T)| > Wc} the (random) time 
when Wc is first exceeded. This time plays the role of the "onset" of large fiuctuations 
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of the piston velocity. One would expect, based on (5.4) that 

T^^hi{WcL/C)l\iiR (5.5) 

i.e. Tc grows as InL when L increases. 

We found Tc experimentally for Wc = 0.1 and Wc = 0.15 and checked that (5.5) 
agreed well with the data, see Fig. 11. By the least squares fit we estimated C = 0.45 
and R — 1.6. 

To summarize our experimental observations, we conclude that the random fluctu- 
ations of the density function pi(y, v,r) and the piston coordinate Y]:^{t) grow expo- 
nentially in time r, and at times r ~ logL they become large even on a macroscopic 
scale. At that point the evolution of the system deviates far from the solution of the 
hydro dynamical equations (H1)-(H4), and they become completely separated afterwards. 

Interestingly, our observations do not indicate that the convergence (1.12)-(1.13) 
claimed in Theorem 1.2 fails on any interval of time. In fact, if the random fluctuations 
behave as CK^/L, as predicted by (5.4), then (1.12)-(1.13) should hold as L ^ oo 
on every flnite interval (0,r*). However, a rigorous proof of this fact would be a very 
challenging task. 

Next, we also examined numerically and heuristically how the system behaves asymp- 
totically, as T — > cxo. On physical grounds [B], wc expect the system to approach thermal 
equilibrium, see Section 1, i.e. the velocity distribution of gas particles should converge 
to a Maxwellian. 

We used the Kolmogorov-Smirnov statistical test to verify the convergence of the 
velocity distribution to a normal law. At any given time r > 0, let 

Fr{u) = : Vi < u}/N 

be the empirical (cumulative) distribution function of particle velocities. For the corre- 
sponding normal distribution function ^{x), we compute 

Dt = sup \F^{u) - ^{u)\ 

— oo<u<oo 

Initially, Dq fa 0.245 for our choice of 7ro{v) in (5.1). If the velocities Vi were independent 
normal random variables, then D^- would be of order 0(1/ ^/N) and the product DrVN 
would have a certain limit distribution, see, e.g. [Lu]. In particular, it is known that 
the probability P{Dt-\/N > 1) ^ 0.2. Based on this, we opted to define the time of 
convergence to equilibrium by 

Teq = inf{r > : Dr^N < 1} (5.6) 

We estimated Teq for various L's and found that Tgq ~ aL^ with some constants a,b > 0. 
By a least squares fit to experimental points we found a = 0.18 and b = 2.47. 

The plot of the product S = DrVN versus r is given on Fig. 12 (for a particular run 
with L = 40). It shows that, after an initial sharp drop over the period < r < 20, the 
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statistic S decreases exponentially in r. Another commonly used (and popular among 
experimentalists) statistic to measure closeness to a normal distribution is 

where M2 and M4 are the second and the fourth sample moments of the empirical veloc- 
ity distribution, respectively. Fig. 12 shows that S' converges to zero in a similar manner 
(for the same run with L — 40) . 




Figure 12: In 5* (thick line) and In 5" (thin line) as functions of r. 



The convergence to a thermal equilibrium can also be justified mathematically. If 
one fixes three integrals of motion - the total energy E and the numbers of particles in 
the left and right compartments A^^ and Nji - the dynamics can be reduced to a billiard 
system in a high-dimensional polyhedron by standard techniques, as we show next. 

Let {xi}, i = 1, . . . , Nr, denote the ^-coordinates of the particles to the right of the 
piston, and {xi}, i = — 1, . . . , — A"l, those to the left of it (ordered arbitrarily). Put 
Xq — X\/M, where X is the coordinate of the piston and M is its mass. Then the 
configuration space of the system (in the coordinates Xi, —N^ < i < Nji) is a polyhedron 
Q C IR^+^ (recaU that N ^ Nl + Nr) defined by inequalities 

< X_Nl, ...,X_i< Xq/^/M <Xi,.. .,XNii < ^ 

It is known that the dynamics of our mechanical system corresponds to the billiard 
dynamics in Q, see [CFS]. That is, the configuration point q E Q moves freely and 
experiences specular refiections at the boundary dQ. The velocity vector 

P = q = {v-Ni^, v-i, V\fM, -yi, . . . , Vf^^ 
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has constant length, since ||p|p = 2E = const. Therefore, the phase space of the billiard 
system is = Q x where Sj^ is the A^-dimensional sphere of radius p = \/2E. 

The billiard system has the Liouville invariant measure fi on Ai, which is the product 
of a uniform measure on the polyhedron Q and a uniform (Lebesgue) measure on the 
sphere , i.e. dp = dqdp. The properties of billiard dynamics depend heavily on the 
curvature of the boundary dQ. In our case Q is a polyhedron, hence its boundary consists 
of fiat sides with zero curvature. A prototype of such systems is billiard in a polygon. It 
is well known that (see, e.g., [C]) 

Fact. For billiards in polygons and polyhedra (and hence, for our mechanical model of a 
piston in the ideal gas) all Lyapunov exponents vanish, and so does the measure-theoretic 
(Kolmogorov-Sinai) entropy. 

Systems with zero Lyapunov exponents and zero entropy arc not regarded as truly 
chaotic, but they still may be ergodic. In fact, billiards in generic polygons are ergodic 
[KMS]. Moreover, for many nonergodic polygons, the phase space is foliated by invariant 
subsurfaces on which the dynamics is ergodic. 

Even though there are no similar results, to our knowledge, for billiards in high- 
dimensional polyhedra, one can expect that they, too, have similar properties. That is, 
they are generically ergodic or become ergodic after trivial reductions. In our case, the 
biUiard in Q is, perhaps, ergodic for typical values of M, or else the phase space is fohated 
by invariant submanifolds on which the dynamics is ergodic, and that those submanifolds 
fill Ai pretty densely. In the latter case, one would hardly distinguish experimentally 
between such a nonergodic system and a truly ergodic one. 

Hence, we can assume that our system is ergodic or very close to ergodic in the above 
sense. Then almost every trajectory eventually behaves according to the invariant mea- 
sure p, independently of the initial state. In particular, for any initial gas density and 
velocity distribution (given by the function 7ro(y, f), see Section 1) the hydrodynamic 
regime for a finite L is only valid on a finite interval of time - eventually the system will 
relax to a thermal equilibrium. We expect in fact that in terms of the "macroscopic" 
variables, say, the one particle distribution function, the system will relax to an effective 
equilibrium, as defined by (5.6) in terms of Tcq, which is much smaller than the exponen- 
tially long time (in L) required for the ergodic theorem. So the real question is how does 
this time depend on L. According to our earlier discussion Tc ~ logL and T^q ~ 

At equilibrium, the distribution of coordinates Xi and velocities Vi are determined 
by the Liouville measure p, which is uniform in the phase space. Physically interesting 
(and only observable) are its marginal measures, i.e. projections, on lower-dimensional 
subspaces. The marginal measures of the velocities are approximately normal for large 
N. 

In particular, each individual velocity Vi converges in distribution to a Maxwellian 
(i.e., normal) random variable with zero mean and variance 2E/N = const. The same 
holds for the "piston" component of the velocity, xq — V\/M, hence the piston velocity 
V will be normally distributed with zero mean and standard deviation const/ -\/M = 
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const/L, as L — > oo. In our case V has standard deviation sJV^/L ^ 0.5 /L. This 
conclusion agrees well with Hollcy's theorem 1.1 and our numerical data. 

The equilibrium distribution of the piston coordinate X is also determined by the 
projection of the uniform measure dq on Q onto the xq axis. Before we do that, let us get 
rid of M in the definition of both Q and xq. A simple change of variable X — xo/\/M 
allows us to redefine Q by 

< X_Nl, ■■■,X_i<X<Xi,..., XNn < L 

Furthermore, rescaling Y — X/L and yi — Xi/L gives a new, simpler, definition of Q: 

< y-N^, ■■■,y-i<y <yi,---,yNR<^ 

This is a variation of the so called Brownian bridge. "Integrating away" the variables yi 
yields the following equilibrium density for Y: 

/(F) = cy^^(i-y)^« (5.7) 

for < y < 1, where c is the normalizing factor that can be computed explicitly: 

c-'^ ['Y^'-il-Yf-dY^ ^^^^"^-y- ^„ (5.8) 
Jo ^ ' (A^L + A^R + 1)! 

Asymptotically, as L ^ oo, we have Nj^ ~ and A^r ~ Assume, for simplicity, 

that Nl^ Nr^ N/2 and denote K = N/2, then 

_ (2A: + 1)! 

Put z — (Y — 0.5)V8K, then the density of z is given asymptotically by 
f{z) = ^fo.5 + -^^''fo.5--^^'' 

C / Z \ 



i^VSK \ 2K^ 
: e" 



1 



2ti 

Hence, Y is asymptotically gaussian with mean 0.5 and variance (4A")~^ = (4L^)~^. 
Therefore, in equilibrium 

1 1 

r-0.5 ^ ^ 

Note that this estimate is independent of the piston mass. 
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Next, by using (5.7) one can easily compute the probability that the piston coordinate 
Y deviates from its mean value 0.5 by a fixed amount d > 0, say d = 0.1: this probability 
is < const ■ e~"^ with some a = a{d) > 0. Therefore, we have observed experimentally 
a very rare event whose probability was exponentially small in A^. If we started with an 
equilibrium state, such an observation would be practically impossible. But we started 
with a state described by the density (5.1), which itself has probabihty less than const • 
Q-bN ^ > 0, with respect to the equilibrium measure So we only observed how one 
highly improbable initial state evolved to another highly improbable state along its (very 
slow) transformation to an equilibrium state. It would be interesting to understand why 
the system "chooses" such a strange evolution to equilibrium, i.e. why starting with a 
state described by a "double-peaked" distribution (5.1) the system behaves as a damped 
harmonic oscillator over an extended time interval, with initial oscillations as large as 
1/10 of the system size. 

We are currently working on this problem and will report results in a separate paper 
[CCL]. Below we outline our program and mention some preliminary findings. For 
simplicity, we assume that 7ro(y,i') = 7i'o(|'f|) and ^(0) = L/2, V{0) — 0. Corollary 2.13 
ensures that the hydrodynamical equations have the trivial solution Y[t) = 0.5, W{t) = 
0, and 7i{y,v,t) = 7ro(t/,f) for all r > 0. Now, since the initial configuration of particles 
is generated randomly from a Poisson process with the density 7ro{y,v), the "actual" 
(empirical) density of the particles, such as the one shown on Fig. 10 at r = 0, does 
not exactly coincide with Tro{y,v). Random fluctuations of the empirical density are 
typically of order 0{1/L). Hence, the "actual" initial distribution of particles can be 
thought of as a small perturbation of the function 7ro(?/,f) = 7ro(|f |) and can be written 
as Trod^l) + £7ri(y, v) with e = 1/L and some (random) function 7:i{y,v) of order one. 

Now, we conjecture that the evolution of the mechanical system closely follows the 
solutions of the hydrodynamical equations (H1)-(H4) with a perturbed initial density 
7ro(|w|) + £7ii{y, v), rather than the stationary solution corresponding to the unperturbed 
density 7ro(|f |). This accounts for significant differences between the behavior of the me- 
chanical system and the stationary solution, if the latter is unstable. In particular, two 
trajectories which are initially £-close (in our case e — 1/L) can deviate from each other 
exponentially fast in time r, and at times r ~ — logs = logL will look completely dif- 
ferent. This would be in agreement with our experimental observations and the estimate 
(5.5) of the time Tc ~ logL of the onset of "instability". 

To test our conjecture, we solved the hydrodynamical equations (H1)-(H4) numer- 
ically starting with a perturbed initial density obtained by adding to (5.1) a function 
£-small in the metric (with e ~ 0.01). We found that the corresponding solution 
resembled strikingly well the evolution of the mechanical system described above. In 
particular, the coordinate and velocity of the piston followed large nearly harmonic oscil- 
lations during the interval 10 < r < 30. The corresponding plots of the piston position 
and velocity along the perturbed solutions of (II1)-(II4) were almost indistinguishable 
from our Figures 6 and 7. Hence, the behavior of the mechanical system can be traced 
to that of the perturbed solutions of the hydrodynamical equations, and the instability 
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of the latter becomes an important issue. 

When wc were finishing the present paper, we received a message from E. Caghoti 
and E. Presutti who (a) proved that the hydrodynamical equations (H1)-(H4) are stable 
when TTodt'l) is monotonically nonincreasing in 1^1, i.e. 7rd(|f |) < 0, and (b) suggested that 
they might be unstable for our class on non-monotone 7ro(|f |). We checked the suggestion 
(b) for our particular density (5.1) and found that it was indeed correct; we proved that 
small perturbations grow exponentially in r. 

Conversely, when we simulated a particle dynamics with a nonincreasing initial density 
TTodvl) the oscillations essentially disappeared (to this end we tried a uniform "fiat" 
function given by TTodi'l) = 1 for \v\ < i^max and a triangular one 7ro(|v|) = 1 — l^l/vmax-)- 
On the other hand, the particle velocity distribution still approached a Maxwellian, albeit 
at a somewhat slower pace. 

Finally, we describe some other open problems related to the piston dynamics. 

1. It is clear that recoUisions of gas particles with the piston have a very "destructive" 
effect on the dynamics in the system. However, we need to distinguish between two types 
of recollisions. 

We say that a recollision of a gas particle with the piston is long if the particle hits a 
wall X = OT X = L between the two consecutive collisions with the piston. Otherwise 
a recollision is said to be short. Long recollisions require some time, as the particle has 
to travel all the way to a wall, bounce off it, and then travel back to the piston before it 
hits it again. Short recollisions can occur in rapid succession. 

We have imposed the velocity cut-off (P4) in order to avoid any recollisions for at 
least some initial period of time (which we called the zero-recoUision interval). More 
precisely, the upper bound Vmax guarantees the absence of long recollisions. Without it, 
we would have to deal with arbitrarily fast particles that dash between the piston and 
the walls many times in any interval (0, r). On the other hand, the lower bound i'min was 
assumed to exclude short recollisions. 

There are good reasons to believe, though, that short recollisions may not be so de- 
structive for the piston dynamics. Indeed, let a particle experience two or more collisions 
with the piston in rapid succession (i.e. without hitting a wall in between). This can 
occur in two cases: (i) the particle's velocity is very close to that of the piston, or (ii) 
the piston's velocity changes very rapidly. The latter should be very unlikely, since the 
deterministic acceleration of the piston is very small, cf. Theorem 2.12c. In case (i), the 
recoUisions should have very httle effect on the velocity of the piston according to the 
rule (1.2), so that they may be safely ignored, as it was done already in some earlier 
studies [H, DGL]. 

We therefore expect that our results can be extended to velocity distributions without 
a cut-off from zero, i.e. allowing fmin = 0. 

2. In our paper, L plays a dual role: it parameterizes the mass of the piston (M ~ L^), 
and it represents the length of the container {0 < x < L). This duality comes from our 
assumption that the container is a cube. 
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However, our model is essentially one-dimensional, and the mass of the piston M and 
the length of the interval < a; < L can be treated as two independent parameters. In 
particular, we can assume that the container is infinitely long in the x direction (i.e., that 
L is infinite), but the mass of the piston is still finite and given by M ~ {this L is the 
size of the container in the y and z directions). In this case there are no recoUisions with 
the piston, as long as its velocity remains small. Hence, our zero-recoUision interval is 
effectively infinite. As a result. Theorem 1.2 can be extended to arbitrarily large times. 
Precisely, for any T > we can prove the convergence in probability: 

P I sup \Yl{t,(jj) - Y{t)\ < Ct In L/L) ^ 1 
\o<T<r / 

and 

P I sup \WLiT,cu) - W{t)\ < CTlnL/L] 1 

\0<T<T / 

as L — > oo, where Ct > is a constant and Y{t) and W{t) — Y{t) are the solutions of 
the hydrodynamical equations described in Section 2. 

3. Along the same lines as above, we can assume that the container is d-dimensional with 
(i > 4. Then the mass of the piston and the density of the particles are proportional to 

L'^~^ rather than L^. 

When d is large, the gas particles are very dense on the x, v plane. This leads to a 
much better control over fluctuations of the particle distribution and the piston trajectory. 
During the zero-recoUision interval, for example, the piston trajectory is L~(''~^)/^-close 
to its deterministic trajectory. This is an easy modification of the results of our Sec- 
tion 3. During the one-recollision interval, the piston trajectory is L~'^^'^~^)/^-close to its 
deterministic trajectory. This can be shown with the methods developed in Section 4 
but requires some extra work. Moreover, the methods and results of that section can be 
extended to the A;-recollision interval (r^, r^+i) for any k > 1. It can be shown that there 
is a djfc > 3 such that for all d> dk we have 

PI sup \Wl{t,uj)-W{t)\<L-''] ^1 

as L —>■ oo, here b > depends on k and d. This extension, however, requires quite 
substantial work, which is beyond the scope of this article. The upshot is that a high 
dimensional piston is more stable than a lower dimensional one. 

It would be interesting to investigate other modifications of our model that lead to 
more stable regimes. For example, let the initial density vro(?/, v) of the gas depend on the 
factor a = el!^ in such a way that 7ro(|/, v) = a~^p{y, v), where p{y, v) is a fixed function. 
Then the particle density grows as a — > 0. This is another way to increase the density 
of the particles, but without changing the dimension. One may expect a better control 
over random fluctuations in this case, too. 
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Appendix 

In this section, we derive various probabilistic estimates on the distribution of gas parti- 
cles and their velocities. The number of particles K = Njj in any domain D on the x, v 
plane at time i = is a Poisson random variable, and we need bounds on its large devi- 
ations. Perhaps, some of our estimates are known in probability theory, but we include 
proofs for the sake of completeness. 

Lemma A.l Let K be a Poisson random variable with parameter A > 0. Then for any 
A> X we have 

P{K >A)< e^-^-^M^/A) 

and for any ^ < A 

P{K <A)< e^-^-^MA/A) 
Proof. The moment generating function of K is 

First, let A > A. Then, obviously, for alH > 

'Pxit) > e^* • P(K > A) 

Hence for all t > 

P{K >A)< e^(^*-i)-^* 
The expression on the right hand side takes minimum at 

t = \n(A/X) > 

This proves the first part of the lemma. 
Now let A<X. Then for alH < 

(pK{t) > e^* • P{K < A) 
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Hence for alH < 

P{K <A)< e^(^*-i)-^* 
The expression on the right hand side takes minimum at 

t = ln{A/\) < 

This proves the second part of the lemma. □ 

Lemma A. 2 Let K be a Poisson random variable with parameter A > 0. For any b > 
there is a c> such that for allO < B < b\/X we have 



P{\K-X\ > B\fX) < 2e" 



CS2 



Proof. A direct apphcation of the previous lemma gives 

P{\K -X\> bVX) < 2e-^'^(«) (A.l) 

where 

9^^) - -2 - 2f^sds ^^-^^ 

and q = B/\f\ in (A.l). By direct inspection one can verify that the function g{q) is a 
positive and strictly monotonically decreasing function on the interval < g < oo. We 
complete the proof by setting c = g{h). □ 

Lemma A. 3 Let Aq > and a > 0. Then for all sufficiently large L > and every 
Poisson random variable with parameters X> Xq we have 

P{\K-X\ > aVX InL) < L-'^i"'"^^ 

where d — a-\/Ao/2. 

Proof. Using (A.l) with B — alnL gives 

P{\K - A| > aVX InL) < 26" ('^i^^)'^^"^"^/^) 
with g{q) defined by (A. 2). Observe that, for large L, 

a In L \ a/Ao In In L 



J alnL 



9 

This complete the proof. □ 

Remark. In most of our applications, A is large, of order X ^ L^ with some 6 > 0. 
Hence, the factor In L is small compared to \/X. 
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Corollciry A. 4 Let Aq > and a > 0. Let K be a Poisson random variable with 
parameter A < Aq. Then for all sufficiently large L > we have 

P{K > alnL) < L-^^^^^^ 

where d — a/2. 

Proof. The case A = Aq easily follows from the previous lemma. Now, if A < Aq, then the 
event X > alnL is even less likely than it is for A = Aq. □ 

Next, we need to study another random variable related to a Poisson process. For 
any domain D on the x, v plane consider the sum of the velocities 

Z = Zd= V 

of the particles in D at time 0. We assume that D C {fmin < < ^max} (the case 
D C {— fmax < V < — Vmin} is Completely symmetric and analogous). By projecting the 
domain D onto the v axis we obtain a Poisson process on the interval 

I — (^minj ^max) 

with density 

niv)^L^f p.Mdx 

JDn{u=v} 

Now the random variable Z can be described as follows. 

Consider a one-dimensional Poisson process with density 7r{v) on the interval /. This 
means that for any subinterval J <Z I the number of points in J, call it Nj, is a Poisson 
random variable with mean 

E{Nj) ^ j^T:{v)dv 

Each realization uj of this process is a finite subset of /. Consider a random variable 

We will call Z an integrated Poisson random variable. 

If we fix a large n> 1 and partition I into small intervals 



7n 



i i + 1 



_n n 

i — 0,1,2, . . ., then we can obviously bound Z by 



E^^A.<Z<E^^A. (A.3) 
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where A^A; is the number of points of the process in the interval Aj. Note that A^A; are 
independent Poisson random variables with parameters 



\ = E{N^^) = [ 7T{v)dv 
Therefore, the moment generating function (pz{t) = E{e*^) of Z is bounded by 



exp 



< E{e^^) < exp 



Taking the limit n — > oo we obtain 

ipz{t) = E{e'^) = exp 



e*^ - DttM dv 



By using (A. 3), it is also easy to find the mean value 



E{Z) = Hz = J^V7r{v)dv 



and the variance 



Var(Z) = (J% = j v'^7r{v) dv 
Note that Z is related to a Poisson random variable K — Nj with parameter 

In particular, we have 



(A.4) 

(A.5) 
(A.6) 



(v) dv 



hence 

We also have 



■Umini^ < Z < Vmav-fr 



(A.7) 
(A.8) 



vli^Xz <o\< vl^^Xz (A.9) 



The random variable Z admits bounds on large deviations similar to the ones we 
found for Poisson random variables: 

Lemma A.5 For any b > there is a c> {determined by b, Vram CLnd Vmax) such that 
for any integrated Poisson random variable Z and all < B < baz we have 



P{\Z-fXz\> Baz) < 2e 



-cB^ 
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Proof. Put, for brevity, ji — jJiz and a — oz- We will show that 

P(Z > [x + Ba)< e-"^'' 



(A.IO) 



(the same bound for P{Z < /i — Ba) is proved similarly, as we did that in Lemma A.l). 
For alH > we have 

Mt) > 6^'^+^'^)* ■P{Z> 11 + Ba) 



hence 



P{Z > i^ + Ba) < exp 



J^ie^" - 1)tt{v) dv-{ii + Ba)t 
We substitute t — Bs/a with s > to be chosen later and expand e^'" into a Taylor series: 

P{Z>^ + Ba) < exp [ / ^-^^ .{v) dv - ^ - B\s 

The first two terms with n — 1 and n — 2 give 

r Bsi 
Ji a 

by (A. 5) and 

{Bsv)"" 



Bsv . , , Bsu, 

TT V dv = 



2(72 



t:{v) dv 



B's 



2„2 



by (A. 6), respectively. 
Therefore, 



P{Z > n + Ba) < exp [- ^^(s - s72 - 

7r(v) civ 



n=3 



m a^ 



Assuming B/a < b and using (A. 6) gives 

{D.s/a)"--'j]r"K{r)dv 



E 

n=3 



nl a^ 



2 

min 



< 



7)^ / 

max j sbvn 



- 1 



Now if s is small enough, then we have 



c := s - — - ^ rgs&w _ i) > 
2 V • ^ ^ 



This proves (A.IO), and hence the lemma. □ 
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Lemma A. 6 Let Aq > 0. For all sufficiently large L > and any integrated Poisson 
random variable Z with Xz > Aq 

P{\Z-iiz\ > (JzlnL) < L-''i°i°-t' 

where d > is a constant determined by Aq, t'min, fmax- 

Proof. Let 

^ = — 

V ■ 
^min 

First, if liiL < baz, then the result easily follows from the previous lemma. 
Now, assume that 

In L>baz (A. 11) 

Using the inequalities (A.8)-(A.9) gives 

4 / i^^Lf 



Xz<^< 



2 



min 



b^v"^ ■ 



and hence 

^;inax(ln ly 



I^Z < t'maxAz < — (A. 12) 



and also 



V 



I^Z < Vraa^Xz < -^CFz (A. 13) 



Multiplying (A. 12) and (A. 13) and taking the square root gives 



liz<^azlnL=-azlnL (A. 14) 

Therefore, since Z is a positive random variable, we have 

P{\Z - fXz\>(TzlnL) = P{Z > fXz + (TzlnL) 

< P{Z>azlnL) 

Moreover, combining (A. 14) with (A. 11) and (A. 9) gives 

- (7z In L = (7z In L > ^— cr^ > v^^Xz 

and also, by (A. 9) 

-az\nL> - Vmin\/Az In L 
Now, since Z < v^g^K by (A. 7), we have 

P{Z>azlnL) < P (^Z > v^^Xz + „ ^min 

< p(K>Xz + ^^lnL 

Now the result follows from Lemma A. 3. □ 
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Corollciry A. 7 Let Aq > 0. Let Z be an integrated Poisson random variable with Xz < 
Aq. Then for all sufficiently large L > we have 

P{Z > InL) < L-'^i"^!"^ 
where d > is determined by Aq, Vmin, 'I'max- 

This immediately follows from Corollary A. 4. 
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